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2 
PREFACIO A LA SEGUNDA EDICIÓN 


En la segunda edición este folleto fue sometido a una modifi- 
cación. El método de iteraciones está expuesto aquí en la base 
del concepto de la aplicación contraída lo que hizo posible relatar 
de ella antes de introducir el concepto de la derivada. Se ha 
ampliado considerablemente la parte del libro dedicada a la reso- 
lución aproximada de los sistemas de ecuaciones. Por fin, todos 
los problemas van acompañados de las indicaciones sobre la reso- 
lución de ellas. 


PREFACIO A'LA PRIMERA EDICIÓN 


El objeto principal de este folleto es exponer varios métodos 
de la resolución aproximada de las ecuaciones. El valor práctico 
de estos métodos es indiscutible. No obstante, se les da poca 
atención, como en las escuelas medias tanto en la escuela superior, 
Por esta razón ocurre con frecuencia que un graduado de la escuela 
superior, en la que él cursó las matemáticas, experimenta dificul- 
tades encontrándose ante la necesidad de resolver una ecuación 
transcendente más sencilla. Con la resolución de las ecuaciones 
tropiezan no sólo los ingenieros, sino también los técnicos, obreros- 
racionalizadores y los representantes de muchas otras profesiones, 
El conocer los métodos de la resolución aproximada de las ecua- 
ciones es también útil para los escolares de clases superiores. 

Puesto que la mayoría de los métodos de la resolución apro- 
ximada de las ecuaciones está ligada con el concepto de la derivada, 
nos vimos obligados a introducir este concepto. Tomamos por base, 
para ello, los evidentes razonamientos geométricos. De este modo, 
para leer la presente obra son suficientes los conocimientos de las 
matemáticas en los límites del noveno prado de la escuela secun- 
daria. 

Para el libro fue utilizado el curso de las conferencias dictadas 
por el autor ante los alumnos de las clases 9 y 10, que asistían 
el circulo matemático escolar adjunto a la Universidad Estatal 
M. V. Lomonosov de Mosaí. 

El contenido de estas conferencias fue aplicado por S. 1, Schwar- 
zburg*, profesor de la escuela secundaria No. 425 de Moscú, para 


* Profesor, doctor en ciencias pedagógicas, miembro 
correspondiente de la Academia de las ciencias pedagógicas de la URSS, 
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el trabajo fuera de clase con los alumnos de novena clase. El autor 
expresa su gratitud a S. I. Schvarzburg por haber elaborado y 
prestado a la disposición del autor los problemas que se resuelven 
por el método de iteraciones y que fueron utilizados al escribir 
este libro. 

La profunda gratitud el autor expresa para V. G. Boltyánski 
cuyas observaciones contribuyeron en un grado considerable a refor- 
zar el carácter práctico de la variante original del manuscrito, 
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sto 
INTRODUCCIÓN 


Al estudiar las matemáticas en el curso escolar, mucho tiempo 
se dedica a la resolución de las ecuaciones y de los sistemas de ecua- 
ciones. Primero se abordan las ecuaciones de primer grado y los 
sistemas de éstas. Luego aparecen las ecuaciones cuadráricas, bicua- 
dráticas e irracionales. Y por fin, un alumno se entera de las ecua- 
ciones exponenciales, logarítmicas y trigonométricas. 

Esta atención para las ecuaciones no puede considerarse casual. 
Se explica por la importancia que tienen las ecuaciones para las 
aplicaciones prácticas de las matemáticas. Sea cual fuese la apli- 
cación, para obtener una respuesta definitiva, se debe resolver, con 
mayor frecuencia, las ecuaciones o los sistemas de ecuaciones. 

En las clases escolares recurren frecuentemente a las ecuaciones 
para resolver problemas fisicos. Examinemos, por ejemplo, el 
problema siguiente: 

Una piedra se deja caer en un pozo. Hállese la profundidad 
del pozo, si se conoce que el sonido provocado por la caída de la 
piedra se ha captado T segundos después de darse la piedra al 
fondo del pozo. 

Si designamos por x la profundidad del pozo, para deter- 
minar x obtenemos la ocuación 


donde v es la velocidad de propagación del sonido en el aire 
2x i 4 x è 
GA el tiempo que dura la caida, y La es el tiempo 


que transcurre hasta el momento de captar el sonido) 


Es una ecuación irracional. Haciendo y= = y, reducimos- 
la a la ecuación cuadrática 


—+ 


que se resuelve según la fórmula conocida. 
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K_—_—_ a 


Las ecuaciones también se aplican para resolver problemas 
geométricos. Por ejemplo, el problema de división del segmento AB 
de longitud } en la razón media y extrema (es decir, en los 
segmentos AC y CB tales que AB: AC = AC:CB) conduce a la 
resolución de la ecuación cuadrática 


*+l—P=0, 


donde con x está designada la longitud del segmento AC. 
A una ecuación más compleja se reduce el problema de divi- 
sión del ángulo «a en tres partes, Esta ecuación tiene la forma 


4x? — 3x — cosa =0, 


donde x = cos $. Las ecuaciones de este tipo, llamadas cúbicas. 


no se estudian en las escuelas, pero en el curso del álgebra 
superior se demuestra que la fórmula para resolver las ecuaciones 
aúbicas también existe (véase más abajo la fórmula (3). 

Sin embargo, en la física surgen frecuentemente los problemas 
que conducen a las ecuaciones más complejas para las cuales no se 
dan fórmulas de resolución en la escuela média, ni aun en la 
Universidad. Tomemos, por ejemplo, una barra de hierro (viga, 
como dicen los técnicos) con extremos empotrados. Sometida a un 
golpe, la barra empezará a efectuar oscilaciones transversales, 
Como se demuestra en la fisica matemática, con el fin de hallar 
la frecuencia de estas oscilaciones se debe resolver la ecuación 


a > 
ete a (1 


donde e = 2,71828... 

Las escuelas no dan ningunas reglas para resolver la ecuación 
de esta índole. No se debe pensar que esto se explica por la 
brevedad del programa escolar respecto a las matemáticas. 
No existe, en general, la fórmula para resolver la ecuación (1) 
en el sentido escolar. Precisemos esta afirmación. 

Suele decirse que para cierta ecuación existe la fórmula de 
resolución, siempre que las raices de la ecuación puedan ser expre- 
sadas en términos de las magnitudes, que figuran en ella, por 
medio de las operaciones aritméticas, la extracción de raices o con 
ayuda de las funciones exponencial, logarítmica, trigonométricas 
y trigonométricas inversas. En este sentido la ecuación cuadrá- 
tica x? + px+q=0 tiene la fórmula de resolución que tiene por 


expresión 


(2) 


Existe también la fórmula de resolución para la ecuación 
cúbica ® 
+ px+g=0. 


Ella tiene la forma 


2 PG: pa 2,2 
sole tala lat (3) 


No obstante, la aplicación práctica de la fórmula (3) está 
obstaculizada con toda una serie de inconvenientes y requiere el 
empleo de los números complejos. 

Existe también la fórmula para resolver las ecuaciones de 
cuarto grado, pero ésta es tan compleja que aqui no la vamos 
a indicar, 

Para las ecuaciones de quinto grado y de grados superiores 
el asunto es peor todavia. Abel **), un matemático noruego, demos- 
tró en el año 1826 que para n> 5 no exist ía fórmula que expre- 
sara la solución de la ecuación algebraica 


aox" + ax +... +44,=0 


por medio de las operaciones aritméticas y la extracción de raices, 


Y A esta fórma se reduce cualquier ecuación cúbica 
ax ax mxa 


mediante la sustitución x + ¿2-= y 
lo 
»s! Sobre la vida y el trabajo creador de Abel véase 
el libro de O. Ore “Notable matemático Niels Hendrik Abel”. Los trabajos 
de Abel fueron desarrollados por el matemático frances E, Galois. De éste 
hay un libro “Favorito de los Dioses” escrito por L. Infeld. 
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Pueden haber fórmulas de resolución* sólo en algunos casos 
particulares de las ecuaciones algebraicas cuyo grado es superior 
a cuarto. 

Si los matemáticos se limitaran al estudio de las ecuaciones 
que admiten una solución exacta, es decir, la solución obtenida a 
base de cierta fórmula, podríamos imaginarnos una conversación de 
un ingeniero con el matemático. 

Ingenlero; Al calcular la construcción, he llegado a una ecuación de 
este tipo (le muestra al matemático la ecuación). Tengo que resolverla 
urgentemente, pues dentro de un mes debo tener terminado el proyecto. 

Matemático: Con mucho gusto quisiera ayudarle, pero para las ecua- 
ciones de este tipo no hay fórmula de resolución. 

Ingeniero: ¿Y no podría Ud. deducir la fórmula? 

Matemático: No vale la pena. Hace tiempo se demostró que para estas 
ecuaciones no existia la fórmula de resolución. ` 

Como resultado de semejante conversación la opinión del 
ingeniero sobre los matemáticos y sus posibilidades se empeora- 
ría, probablemente, en el grado considerable. Afortunadamente, 
tal conversación no puede tener lugar. El hecho es que, habitual- 
mente, al ingeniero no le es necesaria la fórmula para resolver 
tal o cual ecuación. Lo que él necesita es la respuesta obtenida 
con el determinado grado de precisión y no le interesa el modo 
por cuyo intermedio se obtiene la respuesta. La propia fórmula 
sólo sirve para que se pueda calcular la respuesta con el grado 
necesario de precisión. 

Supóngase, por ejemplo, que la fórmula está obtenida, pero la 
respuesta, hallada según esta fórmula, se representa en la forma 
x=3+ yi3. Está claro que esta respuesta no puede ser usada 
inmediatamente (sería absurdo pedir que un ajustador haga una 
pieza de longitud igual a 3+ Vi3cm). Para los fines prácticos se 
debe expresar va en el sistema decimal de cálculo, tomando tal 
número de signos decimales que se requiera para el problema 
práctico dado. 

Asi pues, el ingeniero será totalmente satisfecho, si el mate- 
mático le indica tal o cual modo de calcular la raíz de la ecua- 


* De las ecuaciones algebraicas pueden leerse los li- 
bros de la serie “Lecciones populares de matemáticas” (cm ruso): 
A. G. Kútrosh, “Ecuaciones algebraicas de grado arbitrario”, Editorial 
Mir, 1976. 

L R. Shafarevich, “Acerca de la resolución de las ecuaciones de grados 
superiores”, Gostejizdat, 1954. 
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ción con el grado necesario de precisión. En las matemáticas 
está elaborada toda una serie de los métodos para el cálculo 
aproximado de las ecuaciones. Esta obra se dedica precisamente a 
la descripción de algunos de estos métodos. 


$2 
APROXIMACIONES 
SUCESIVAS 


La mayor parte de los métodos de la resolución aproximada 
de las ecuaciones se basa en la idea de las aproximaciones suce- 
sivas. Esta idea se usa para solucionar no sólo las ecuaciones, 
sino también algunos problemas prácticos. 

El método de aproximaciones sucesivas se emplea, por ejemplo, 
por los artilleros. Si se necesita batir algún objetivo, en la pieza 
de artillería se establecen las graduaciones correspondientes de ele- 
vación- y deriva, después de lo cual se efectúa un disparo. Fallado el 
tiro, a base de la observación del punto de explosión del proyectil 
se introducen correcciones en las graduaciones de elevación y deriva, 
después de lo que se efectúa el disparo siguiente. Realizadas unas 
cuantas aproximaciones, la elevación y la deriva se gradúan de 
tal manera que el objetivo resulta derribado. 


Ao 4% 
As 


Fig. 1 


A veces las aproximaciones sucesivas se necesitan para determi- 
nar el punto de puntería, Supongamos que un cañón antiaéreo, 
ubicado en el punto O, efectúa el tiro contra un avión en vuelo 
(fig. 1). Si apuntamos el cañón al punto Ay, donde se encuentra 
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el avión en el instante dado, tendremos tiro fallado: mientras el 
proyectil vuela, el avión se desplazará al otro punto A,. Este punto 
se determina de una manera bastante fácil, si se conocen las velo- 
cidades del proyectil y del avión. No obstante, si dirigimos el 
proyectil al punto A,, de nuevo podemos errar el blanco. En efecto, 
como resultado del cambio en la inclinación del tubo del cañón 
se altera la ley de movimiento del proyectil. Por esta razón, el 
tiempo que tarda el proyectil en cubrir la distancia OA, se diferen- 
cia del tiempo que se necesita para alcanzar el punto A,. De 
resultas, el proyectil no dará en el blanco. Sin embargo, el error 
será menor en comparar con el caso de puntería en el punto Ap. 
Para que el error sea aún menor, hace falta determinar el tiempo 
durante el cual el proyectil cubre la distancia OA, y, además, 
hallar el punto en que, al pasar este lapso, se ubicará el avión. 
Este punto A, será la aproximación siguiente para el punto de 
puntería buscado. Después de esto se debe determinar el tiempo 
durante el cual el proyectil alcanza el punto Az y calcular el punto 
Aj en el que se ubicará el avión, pasado este tiempo. Realizadas 
unas cuantas aproximaciones, encontraremos, con el grado necesario 
de precisión, el punto de punteria. 

El método de aproximaciones sucesivas se emplea también para 
la resolución de muchos otros problemas. 

Supongamos, por ejemplo, que se necesita transportar la arena 
desde las canteras de arena A;, ..., A, hacia varias obras de 
construcción B,, ..., Bw. Sea a, el rendimiento de la cantera A, 
medido en toneladas por día. Supongamos que la obra B, consume 
b, toneladas por día y que el coste de la transportación de una 
tonelada de arena desde la cantera A, hacia la obra B, es igual 
a Ca (esta magnitud depende de la distancia entre A, y Bj, 
el estado de caminos, etc.). 

Con el fin de elaborar el plan del transporte, compongamos 
la tabla 1. En esta tabla con xy está designada la cantidad de 
arena que se traslada desde la cantera 4, hacia la obra B,. 


Tabla 1 
B, T B, E Bm 
A Xu Ma Xim 
A; Xu Xa Xim 
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Está claro que los números xy deben satisfacer las siguientes 
correlaciones 
Xy +xXa+--«+Xim<4) 


(desde la cantera A, no puede transportar más que a, toneladas 
de arena por día), 
Xir Xa t.e aa = Os 


(a la obra*B, se debe llevar b, toneladas de arena por día). 
Según el plan proporcionado por la tabla 1, el coste del 
transporte de arena es igual a 


C=CjX1 + Cr2X12 +++ H CinXin + 
+ Ca1X21 + C22X22 +-+- + Cinan + (4) 


+ CmiXmi + Cm2Xm2 + --+ E Cae 


El plan se debe elaborar de una manera tal que el valor de 
C sea mínimo. Con este fin se compone primeramente un plan 
preliminar. Se procede, por ejemplo, de tal modo: 

A la cantera de arena A, se le atribuye una obra de cons- 
trucción más próxima a la cantera. Si esta cantera produce más 
arena que se necesita para la obra dada, se le agrega a la cantera 
otra obra, que entre las restantes obras es más próxima a la misma 
cantera. Realizados unos cuantos pasos, el rendimiento de la cantera 
A, resultará agotado. Después de esto, a la cantera A, se le concede 
aquella de las obras restantes que es más próxima a Az, etc. 
Al fin y al cabo, cada obra resultará atribuida a alguna de las 
canteras. 

Sin embargo, el plan elaborado de este modo no es el mejor. 
puesto que al final quedarán unas pocas obras y ellas pueden 
estar muy alejadas de las canteras restantes. Por eso nos veremos 
obligados a cambiar el plan elaborado. Algunas de las obras designa- 
das a las canteras de pequeños números serán tomadas de éstas 
últimas y atribuidas a las otras canteras de números mayores. 

Los métodos de la alteración del plan, que tienen por objeto 
disminuir el coste del transporte, se consideran en el apartado de 
las matemáticas, llamado programación lineal *. 


~ Sobre la programación lineal véase “Introducción 
en la Algebra lineal y la programación lineal” por A. S. Solodovnikov. 
Editorial “Prosvescheniye”, Moscú, 1966. 
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Después de unas cuantas aproximaciones, realizadas en con- 
formidad con estos métodos, obtendremos un plan según el cual 
la suma (4) será minima o poco se diferenciará de la mínima. 

En general, al elaborar un plan (o bien un horario) cual- 
quiera, se toma al principio una aproximación basta la que, 
después, se mejora de manera sucesiva logrando, al fin y al cabo, 
el resultado requerido, 

El maquinado de una pieza en los talleres de una planta 
puede considerarse también como la aproximación sucesiva a la 
forma requerida. Primero se toma una aproximación basta, esto es, 
una pieza de fundición o cualquier otra pieza bruta. Esta pieza es 
trabajada en un torno y se le comunica una forma próxima a 
la deseada. Luego, la pieza se transmite al otro torno de mayor 
precisión, Realizadas unas cuantas etapas del maquinado, es decir, 
Unas cuantas aproximaciones, aparece, por fin, la pieza que deseá- 
bamos obtener. 

La idea general del método de aproximaciones sucesivas está 
bien expresada en uno de los versos de Piet Hein, un sabio y poeta 
danés, los cuales él mismo denominó “gruks”: 


“Espinoso es el camino a la verdad, 
Pero, siguelo y no te demores, 
Fallo, fallo, otra vez failo, 

Pero cada vez menos y menos”. 


$3 
AQUILES Y LA TORTUGA 


Por primera vez las aproximaciones sucesivas se mencionan 
en los razonamientos del Zenón de Elea, filósofo griego, que vivió 
500 años antes de nuestra era. Este filósofo trataba de demostrar 
la ausencia del movimiento en la naturaleza (recuérdense los versos 
de Pushkin llamados “Movimiento”: “El movimiento no existe, 
dijo un sabio barbudo..”). La demostración suya fue como sigue: 
si Aquiles, un hombre más rápido de todos los griegos, hubiera 
tratado de alcanzar corriendo a una tortuga, nunca podría hacerlo. 
Efectivamente, sea la distancia entre Aquiles y la tortuga 1000 
pasos y supongamos que Aquiles hace corriendo 10 pasos por 1 se- 
gundo, mientras que la tortuga cubre 1 paso por segundo. Trans- 
curridos 100 segundos, Aquiles habrá recorrido 1000 pasos que le 
distan de la tortuga, pero durante este tiempo la tortuga se alejará 
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a 100 pasos. Transcurridos 10 segundos más, Aquiles cubrirá co- 
rriendo estos 100 pasos, pero la tortuga estará 10 pasos avanzada. 
Para superar estos 10 pasos Aquiles necesita 1 segundo, durante 
el cual la tortuga se desplazará un paso más. De este modo, 
razonada Zenón, la tortuga siempre irá delante y Aquiles nunca 
la alcanzará. Por consiguiente, el movimiento no existe, 

Este razonamiento de Zenón es, por supuesto, nada más que una 
paradoja graciosa. El movimiento es una propiedad imprescriptible 
de la materia. Es que los versos de Pushkin sobre Zenón son 
como sigue: 


“El otro viejo se puso a andar silencioso, 
No podría ser la objeción más curiosa, 
Elogiaban todos la respuesta ingeniosa”. 


Todo escolar calculará sin dificultad el tiempo que es necesario 
para que Aquiles alcance la tortuga. Para ello seria suficiente 
escribir la ecuación 


10x — x = 1000, i (5) 


en la que con x está designado el tiempo buscado. 
De esta ecuación se obtiene 


1000 Le 
x= s= lllğs 


No obstante, el razonamiento de Zenón puede considerarse 
como un método originai de la resolución aproximada de la 
ecuación (5). 

Para concretar, en la ecuación (5) pasemos x al segundo miembro 
y dividamos ambos miembros de la ecuación entre 10. Obtendre- 
mos la ecuación 


=10+ i (6) 


Si prescindimos en el segundo miembro del sumando 5 (es 


pequeño al comparar con x), obtendremos para x un valor apro- 
ximado x, = 100. Ahora podemos precisar la respuesta, sustituyendo 
en el segundo miembro la aproximación hallada de x, = 100 en 
lugar de x. Obtendremos para x el valor más exacto, a saber, 
xz = 100 + 10 = 110. Sustituyendo este valor en el segundo miembro 
de la ecuación, encontraremos la aproximación siguiente x, = 100 + 


2-7 
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+ ¿lO 111. De este modo obtenemos las aproximaciones 


10 
xı = 100, x; = 110, x3 = 111, x4 =111,1 ... 


es decir, los mismos mimeros que se obtenian en el razonamiento 
de Zenón. Estos números están ligados entre si por la correlación 
Xa 


Sata = 100 + E, 0 
que permite calcularlos uno tras otro. Con el aumento de n ellos 


à m 1 i 
se aproximan a la solución exacta, x = Hig» de la ecuación (5). 


El método que acabamos de exponer ha sido coronado por 
re x h 
éxito, porque el sumando Jg *2 Pequeño comparando con x. 
De lo contrario, obtendriamos los números que no se aproxi- 
marían a la solución buscada. Supongamos, por ejemplo, que 
Aquiles compite no con una tortuga lenta, sino con un antilope 
que es capáz de correr con la velocidad de 20 pasos por segundo. 
Con el fin de determinar el tiempo necesario para que Aquiles 
alcance al animal, hace falta resolver la ecuación 


10x — 20x = 1000. (8) 


La solución es x= — 100. Esto significa que Aquiles y el anti- 
lope estaban al lado hace 100 segundos y en el momento dado 
el animal está a la cabeza de Aquiles y la distancia entre ellos 
sólo va a aumentar en lo ulterior. 

Trataremos de resolver la ecuación (8) del mismo modo que la 
ecuación (5). Pasemos el sumando 20x al segundo miembro y 
dividamos los dos miembros de la ecuación entre 10. Obtendremos 
la ecuación 

x= 100 + 2x. (9) 


Hagamos en el segundo miembro xo = 0. Resulta que x, = 100. 
Sustituyendo este valor en el segundo miembro de la ecuación (9), 
obtenemos la aproximación siguiente x, = 300. Continuando este 
proceso, obtendremos los números 


Xo =0; xı =100; x, =300; x3=700, ... 


Vemos que estos números no tienden a la solución exacta, x; = — 100, 
de la ecuación (8). 
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5. 
DIVISION 
EN LAS COMPUTADORAS 


Es posible que el lector de este libro quede asombrado: 
¿cuál es la razón de resolver la ecuación (5) por el método de 
aproximaciones sucesivas, si se puede obtener la respuesta exacta 
de manera muy sencilla? Nos interesa, por supuesto, no la propia 
ecuación (5), sino el método de aproximaciones sucesivas el cual 
nos proponemos aplicar en adelante a las ecuaciones más complejas. 

Por lo demás, hay que decir que el empleo del método de 
aproximaciones sucesivas para resolver las ecuaciones parecidas a 
la (5) se necesitó hace muy poco tiempo, cuando aparecieron las 
computadoras de acción rápida. Existen modelos de las máquinas 
que son capaces de ejecutar sólo tres operaciones aritméticas: 
adición, sustracción y multiplicación. Además, ellos permiten dividir 
entre los números del tipo 2”. ¿De qué manera estas máquinas 
hacen la división entre números cualesquiera? 

Dividir el número b entre el número a significa resolver la 
ecuación ax = b. Puesto que la máquina puede multiplicar y dividir 
entre 2", se puede considerar que *'/, <a < 1 (de lo contrario mul- 
tipliquemos o dividamos los dos miembros de la ecuación ax = b 
entre la potencia correspondiente del número 2). 

Escribamos la ecuación ax =b en la forma 


x=(1-a)x+b. (10) 
Tomemos por la primera aproximación del número x el valor 
xı = b. Designemos con a, el error de esta aproximación, es decir, 


supongamos que x; + a, = b/a. Entonces, de la ecuación (10) obte- 
nemos 


a+ =(1-al+a)+o=(1ax +o+ (aa. (11) 


Como */, <a < 1, resulta que 
0<1l-a<!,. 
Puesto que el coeficiente 1 — a es comparativamente pequeño, en el 
segundo miembro de la ecuación (11) despreciemos el sumando 


(1 — aja,, el cual es, por lo menos, dos veces menor que 4. 
Obtendremos: 


+0 =(1-a)x, +b. 


2 
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El número 
x=(l-0x+b 
lo tomaremos por la siguiente aproximación para x. 


Designemos por az el error de la aproximación xz, es decir, 
hagamos x, + a = b/a. En este caso, de la ecuación (10) se infiere 


que x +% =(1-—d)x2+b+(1— aja. ` 
Despreciando en el segundo miembro de esta ecuación el 
sumando (1 — a)æz, obtendremos la igualdad aproximada 
Xa + 0% 3 (1 —a)x, + b. 
Así pues, a título de la siguiente aproximación se puede tomar 
x = (1 —a)xz + b. 


Siguiendo con este mismo procedimiento, encontraremos que la 
siguiente aproximación tiene la forma 


Xa = (1 — a) x, + b, 


etc. Los números X4, Xz, ..., Xm +...» calculados sucesivamente 
según la fórmula 


Xavi = (1) x, + b, (12) 


van aproximándose al número b/a. Más, en esta fórmula se 
utilizan sólo las operaciones de adición, sustracción y multiplica- 
ción lo que significa que la máquina puede operar conforme a ella. 

El método de división descrito se basa, de hecho, en la fór- 
mula para la suma de una progresión geométrica infinita decre- 
ciente. Para concretar, escribimos la fracción b/a en la forma 


b_ b 
a 1-0-4 


Pero, según la fórmula citada tenemos: 


b 
Ia- g7ttbl -a+b +. 


+b(l- a)" +... (13 


Designemos con x, la suma de n primeros términos de esta 
progresión: 
x =b+b(1-—a)+...+b{1 - a)" ! 
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Es evidente que 
Xas =b+D(1 a) +... + b(1-a)"= 
=b+ (1-a) [b+b(1-a) +... +b -~ a)"™']= 
=b + (1 — a) Xn- 


Esta fórmula y la fórmula (12) coinciden. De este modo, sus- 
tituyendo la fracción b/a por el valor aproximado de x,, sustituimos 
la suma infinita en la fórmula (13) por la suma de los primeros 
n términos. Al aumentar el número de los sumandos n, esta suma 
va aproximándose a la suma de toda la progresión (el decreci- 
miento de la progresión (13) se deduce de lo que '<a<l, 
y, por tanto, 0 < 1 — a < '/3) 


$5 
EXTRACCIÓN DE RAÍCES 
CUADRADAS POR EL MÉTODO 
DE APROXIMACIONES SUCESIVAS 


Mostremos ahora, cómo se emplea el método de aproxima- 
ciones sucesivas para la extracción de raíces cuadradas. En las 
escuelas se estudia un método de extracción de raíces cuadradas 
que permite hallar, uno tras otro, los signos de la raíz buscada. 
Se lo puede considerar también como el método de aproxima- 
ción sucesiva hacia la respuesta. Sin embargo, es bastante compli- 
cado y muy frecuentemente los alumnos lo usan mecánicamente, 
sin darse cuenta completa de lo esencial en la operación. Vamos 
a describir aquí el otro método de extracción de raíces cuadradas 
que se utilizaba ya en Babilonia" antigua mucho tiempo antes de la 
era nuestra. Lo empleaba también Herón, un matemático de Alejan- 
diría. Después, este método fue abandonado, pero en nuestros tiempos 
se le recurre a veces para extraer raices cuadradas en las máquinas 
de cómputo. 

Supongamos, por ejemplo, que se necesita extraer una raíz 
cuadrada del número 28. Elijamos primero algún valor aproximado 
de esta raíz, haciendo, por ejemplo, x, =5. Designemos por a, 
el error de este valor aproximado, es decir, hagamos 
YB=5+%. Pára hallar el valor de «,, elevemos al cuadrado 
los dos miembros de esta igualdad. Obtendremos: 


28 =25+ 100, +a}, 
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Siddi a? + 10a, —3=0. 

Obtuvimos, de este modo, una ecuación cuadrática para a. 
Si tratamos de resolver esta ecuación con exactitud, obtendremos 
el valor de &; = -5 + va. Asi pues, para hallar el valor exacto 
dé a, se debe calcular v28. Pareciera que estamos dentro de un 
círculo vicioso: para hallar 128 se necesita calcular %,, y para 
hallar æ, se debe calcular vos. 

La siguiente consideración viene a ayudarnos. El error a, del 
valor aproximado x, =$ no es significante; de antemano se sabe 
Que es menor que 1. Es menor aun el número a?, Por esto, 
hagamos una tentativa de hallar el valor aproximado de Q,, 
despreciando en la igualdad (14) el sumando pequeño aí. En este 
caso, para a, se obtiene la ecuación aproximada 10a; — 3 x 0. 
De ella se deduce que a, = 0,3. 

Por consiguiente, el valor aproximado de la corrección a, lo 
tenemos. Puesto que V8 =5 4+ a, la segunda aproximación xz 
para |/28 tiene por expresión 

Xx = 5 +0,3 = 5,3, 


Repitamos el procedimiento descrito con el fin de hallar la 
aproximación para v2 más exacta todavía. A saber, designemos 
con az el error del valor xz=5,3, es decir, hagamos 
VB =x, + 0%. Elevemos ambos miembros de esta igualdad al 
cuadrado y despreciemos el pequeño sumando «3. Obtendremos: 
28 = xå +2x20%,, y, por lo tanto, 

28 = x?, 

DER > 

Esto significa que la tercera aproximación para v28 se expre- 
sa mediante la fórmula 


az 


28- x2 _ 28+x} 

Dep w 
Puesto que x, = 5,3, de aquí se infiere que x, = 5,2915.. De un 
modo igual, partiendo del valor aproximado x, = 5,2915, encontra- 
mos la siguiente aproximación x4, expresada mediante la fórmula 


28 + x$ 
= i AS. 


X3 = X3 + 


X4 
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En general, una vez hallado el valor aproximado x, para 
y28, la aproximación siguiente tendrá por expresión 
28 +x? 
y RE ox 
Con ello, todo paso consecutivo del procedimiento conduce a las 
aproximaciones para v28, que son cada vez más exactas. El proceso 
de las aproximaciones se da por terminado cuando la diferencia 
entre Xa+ı Y Xa se hace menos que la precisión predeterminada 
de los cálculos. Por ejemplo, si es preciso calcular /28 con un 
error inferior a 0,0001, será suficiente tomar cuatro aproximaciones 
y hacer V8= 5,2915 (en efecto, x, = 5,2915.. Y X= 5,2915... 

De la manera sumamente igual se extrae la raíz cuadrada de 
cualquier número positivo. A sáber, calculando Ya, eligimos algu- 
na aproximación inicial x,, y luego determinamos las aproxima- 
ciones siguientes, valiéndonos de la fórmula 
a+xi 
2" 

La fórmula (16) se puede deducirla también de los razona- 
mientos algo diferentes de los que hemos aplicado al extraer 
vz. Supongamos que ya hayamos hallado la n-ésima aproxima- 


ción x, para //a. Dado que yg = k, entonces |/a es el 
q CR 
n 


promedio geométrico de los números x, y z. En calidad del va- 


Xa+ = 


(16) 


Xp+1 = 


n 
lor aproximado para este promedio geométrico tomemos el pro- 


e ue ` a 
medio aritmético de los números X, y En es decir, hagamos 
n 


1 a x +a 
Xar 3 At = » 
n 


Esta es precisamente la fórmula (16). 

De este modo, el método aproximado de la extracción de rai- 
ces cuadradas escrito más arriba consiste en lo que en cada paso 
de las aproximaciones el promedio geométrico de los números X„ y 


= se sustituye por el promedio aritmético de los mismos. 


Aclaremos ahora: ¿es que siempre el proceso de las aproxima- 
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ciones sucesivas leva al objetivo al extraer raíces cuadradas, es 
decir, siempre será el asunto tal como lo observábamos cuando. 
Aquiles corría tras la tortuga o, quizás, a veces el asunto recordará 
su competición con el antilope? En el primer caso los matemá: 
ticos dicen que el proceso de las aproximaciones es convergente, 
en el segundo caso el proceso diverge. Vamos a demostrar que al 
extraer las raíces cuadradas nunca nos tropezamos con compli- 
caciones algunas: el proceso de las aproximaciones siempre converge, 
siempre lleva al objetivo. pa 

Para ello comparemos los errores 0, =]/4—xX, Y Uni = 
=|/a — xp+, de dos aproximaciones que siguen una tras otra. 


Valiéndonos de la fórmula (16), el error „+, puede ser escrito 
en la forma 


~- xit+a xi — 2x, a + a 
zarji A 
z 3 
Pero, 
x? -2x, Ya + a = (x, — 0 = o? 
y por esto 
2 
A > an 
r 


Consideramos sólo las aproximaciones positivas x, para Va. 
Por esta razón, de la igualdad (17) se puede sacar la deducción 
de que todos los errores %,, %y, ..., Om ».. Son negativos. En otras 
palabras, todas las aproximaciones Xx», a partir de la segunda, 
son aproximaciones por exceso *; la primera aproximación x, puede 
ser tanto por exceso como por defecto. 

Con ayuda de la fórmula (17) se demuestra con facilidad que 
el valor absoluto del error de la aproximación aproximado de Xx, 
disminuye con cada paso al menos en dos veces. En efecto, la 
igualdad (17) se puede escribir así: 


r 
Xa- Va f 
A (3 Ee) 


ZE 


Por esto A 
Va 


Aia 


lanl- (18) 


dansal 


+! Esto se explica por el hecho de que el promedio 
aritmético es siempre mayor que el promedio geométrico 
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Pero, como x, > 0, entonces 


e M 
2 HA 2 


Por otra parte, según lo mostrado más arriba, para n>2 
tenemos x, >|/a, y por tanto 


1 a 
na” 


De aquí se deduce la desigualdad 


3 A (9) 
2 2% 2 
Explorando las correlaciones (18) y (19), nos convencemos de que 
1 
lanti] < 3il. 


Esto precisamente demuestra la validez de nuestra afirmación: 
con cada paso de las aproximaciones el valor absoluto del error 
va disminuyendo por lo menos en doble, De esto se deduce que 
después del segundo paso de las aproximaciones el valor absoluto del 
error disminuirá al menos en cuatro veces, después del tercer 
paso, al menos en ocho veces, etc. Está claro que al crecer n, 
el valor absoluto del error a, = Va = x, va disminuyendo y tiende 
a cero. Pero, esto significa que los números x, tienden a 

a cuando n crece. 

Ahora pondremos en claro, cómo influye la elección de la 
aproximación inicial x, en el proceso de las aproximaciones, En pri- 
mer lugar indiquemos que esta elección no tiene influencia alguna 
en el resultado final. Es que ya hemos demostrado que con cual- 
quier aproximación inicial x, los errores da, ..., Op de las 
aproximaciones ulteriores tienden a gero cuando n— œœ. De este 
modo, si la precisión necesaria de cálculos está prefijada, se ob- 
tendrá el mismo valor para va con la precisión predeterminada, 
cualquiera que sea la aproximación inicial x,. Si incluso la elec- 
ción de la aproximación inicial no es acertada, llegaremos por fin 
al resultado correcto. Realizados 10 pasos de la aproximación, 
el valor absoluto del error disminuirá más que en mil veces 
(210 = 1024 x= 1000); realizados 40 pasos, por lo menos en mil mi- 
llones (10°?) veces. Por ejemplo, si, al calcular va ponemos 
xı = 106, entonces a, = 10%, y por tanto, «so |< 1076. En otras 
palabras, al principio del proceso el error era cerca de un mi- 
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llón; al final, su valor absoluto se hizo menos que una milio- 
nésima. 

No obstante, la elección de la aproximación inicial influye en 
la duración del proceso de aproximaciones. Si la elección de la 
aproximación inicial no es acertada, tendremos que esperar mucho 
hasta que la diferencia entre x,+; y Xx, se haga menos que la 
precisión dada de los cálculos. La elección acertada de la apro- 
ximación inicial contribuye considerablemente a la aceleración del 
proceso. Por esta razón, se procede con frecuencia de la manera 
siguiente: la aproximación inicial x, se toma de las tablas de 
raices cuadradas y la fórmula 

e A á (20) 
EA 
sólo se emplea para precisar el valor obtenido. 

Este procedimiento es de comodidad particular porque la 
velocidad de disminución del error crece considerablemente cuando 
x, va aproximándose a )/a. En efecto, al deducir la desigualdad 


1 
loas |< llo 
1_ya 
al por el 


número 3 Pero, si x, es próximo a Va, la fracción 


1 ya 


en la fórmula (18) hemos sustituido el factor 


ja 
tps será muy pequeña, por lo que |%,+1|= 
1 2 


es considerablemente inferior a zlat 

Precisemos esta afirmación. Examinemos, a la par con el error 
absoluto EN) = | ya- xl, un error relativo f, del valor aproxi- 
mado x,, es decir, la razón del error absoluto |a,| al valor exacto 
de la raíz Va. Este error tiene por expresión 


lani Xa 
„= — =|- |. 
2 ya ya 


De la igualdad (17) obtenemos la siguiente fórmula para la 
magnitud P,+/: 
ls] _ lanl? 


ya ya 


Bari = 


Puesto que X, > Va, tenemos: 


Así pues, para los errores relativos B, de los valores apro- 
ximados se cumple la desigualdad 


2 
Basi < E, (Q1) 


Por ejemplo, si el error relativo de la aproximación x, es 
igual a 0,01, entonces para la aproximación Xp+, él no supera 
a 0,00005, y para X„+2 no es superior a 0,00000000013. Vemos que 
la precisión de las aproximaciones crece cada_vez más rápido. 
Se puede mostrar que estando ya cercanos a Va, con cada apro- 
ximación ulterior el número de signos exactos se dobla. 

Ejemplo. Calcúlese y238 con un error inferior a 0,00001. 

Rigiéndonos por las tablas de 1. N. Bronstein y K. A. Semen- 
diaev*, encontramos y238 = 15,43. Tomemos 15,43 por x, y ha- 
llaremos x, según la fórmula 
15,43? + 238 

30,86 


Evaluemos la precisión de la respuesta obtenida. Como el error del 
valor 15,43 no supera a 0,01, entonces œ, =0,01, y por eso 
0,01 


PB = 154 < 0,001. 


= 15,42725... 


x= 


Pero, en este caso 


2 
p2 < ooi = 0,0000005. 


Esto significa que el error absoluto de la aproximación x, 
no es superior a 15,43-0,0000005 < 0,00001. En otras palabras, 
en el valor 238 = 15,42725... son exactos todos los siete signos. 

Si deseáramos obtener catorce signos exactos, la tercera apro- 
ximación ya nos daría la respuesta requerida. Sin embargo, esta 
precisión es casi siempre supérflua. 

En conclusión señalemos la siguiente particularidad del método 


* I, N. Bronstein y K. A. Semendiacv, “Manual de 
Matemáticas para ingenieros y estudiantes”, Editorial Mir, 1977. 
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de aproximaciones sucesivas. Al emplear el método ordinario de la 
extracción de las raices cuadradas, cualquier error cometido en 
alguna operación hace completamente inválidos los cálculos ulte- 
riores. Otra cosa es, cuando se usa el método de aproximaciones 
sucesivas. Supongamos, por ejemplo, que como resultado del error 
cometido, en lugar del valor correcto de la n-ésima aproximación 
Xx, Obtuvimos el valor falso y,. Entonces, todos los cálculos ulte- 
riores pueden considerarse como obtención de Va con la apro- 
ximación inicial y, Ya hemos visto anteriormente que cualquiera 
que sea la aproximación inicial, el método de aproximaciones 
sucesivas lleva, al fin y al cabo, al valor de |a con el grado 
deseado de precisión. De este modo, el error cometido va a tender 
en lo ulterior a cero. El único inconveniente que afecta la opera- 
ción consiste en lo que nos veremos obligados a tomar unas cuantas 
aproximaciones de sobra. 

Gracias a la particularidad descrita del proceso de aproxima- 
ciones, al principio del procedimiento se puede calcular las aproxi- 
maciones con menor precisión y tomar la precisión dada sólo para 
las últimas aproximaciones. Esta circunstancia reduce el tiempo 
necesario para los cálculos. 


$6 
APLICACIÓN DEL MÉTODO 
DE APROXIMACIONES SUCESIVAS 
A LA EXTRACCIÓN DE RAÍCES CON 
EXPONENTE NATURAL 


El método de la extracción de raíces cuadradas expuesto más 


arriba, se lo puede emplear también para extraer raíces con cual- 
quier exponente natural. Para ello sirve la fórmula siguiente” 


(xata ktat., (22) 
en la que con los puntos están designados los términos que con- 


tienen a, al, etc. 
Demostremos esta fórmula. Del curso escolar sabemos que 


(x +a} =x? 2x0, 
(x+ aj? = x? data dra ad, 


* Esta fórmula se deduce del binomio de Newton, 
pero no se supone aquí dar a conocer al lector la fórmula del binomio. 
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Se puede escribir estas igualdades en la forma siguiente 
(x +a? =x 2x0 ++... (23) 
(x+? e Ha.. (24) 


De este modo, la fórmula (22) está demostrada para k = 2 y k = 3. 
Multipliquemos, ahora, ambos miembros de la fórmula (24) por 
x +a. Obtendremos: 

{x +a) = (x? + 3xa +.. xt 


Si en esta expresión abrimos paréntesis, se obtendrán un su- 
mando x* que no contiene a y dos sumandos, 3x°x y x'a, 
en los que a interviene en la primera potencia, en cuanto a los 
sumandos restantes, éstos contienen «a en las potencias segunda y 
superiores. Por eso, se puede escribir que 


(x+ at= xt 3a a = a (25) 


(igual que antes! con los puntos están designados los términos que 
contienen a?, a°, etc.). 

Asi pues, la fórmula (22) queda también demostrada para k = 4. 
Del mismo modo, de la fórmula (25) se deduce que 


Garat to... (26) 


Es evidente que siguiendo este procedimiento se puede demos- 
trar la fórmula (22) para cualquier exponente k entero y positivo. 
Volvamos ahora a la extracción de la raiz de cualquier expo- 
nente natural k. Supongamos que ya está encontrada alguna apro- 


di r Sad 4 
ximación x, para la raíz buscada |/a. Designemos con a, el error 


ho 
de esta aproximación, es decir, supongamos que x, +&; = |a. 
Entonces, (x, + )* = a. Pero, valiéndose de la fórmula (22), se puede 
escribir esta igualdad en la forma siguiente: 


x + kx a+. .=0, 


donde con puntos están designados los términos que contienen 
ai, aj, etc. 
1 Aj, ete. 
Si la aproximación elegida x, es suficientemente cercana a 
k 
Va, entonces el error x, de esta aproximación es pequeño, por lo 
que se puede despreciar los términos que contienen las potencias 
superiores del error. Obtenemos, de esta manera. una igualdad 
aproximada 
A+kla a, 
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De esta igualdad se infiere que 
a-% 


A E, 
EA 
£ 
y, por ello, a título de la aproximación ulterior para Va se puede 
tomar el número 
ax a+(k-Dxk 


X= x + = 
A kx! 


Del mismo modo, partiendo de la aproximación x, hallamos 
la aproximación siguiente: 


de a+(k—1)x% 


kx! 
Y, en general, una vez hallada la aproximación x, para 


k 4 LA A . . 
Va, la aproximación posterior se definirá por la fórmula 


Ssa HADS (27) 


Igual que en el caso de la extracción de raíces cuadradas, 
se puede mostrar que el proceso citado converge con elección 
cualquiera de la aproximación inicial x, (con tal de que esta apro- 
ximación sea un número positivo). En otras palabras, cualquiera 
que sea la elección de x,, los números Xi, Xz, ..., Xp... Se 


tr . ; 
aproximan a |'a. El proceso de aproximaciones se lleva a cabo 
hasta que coincidan, dentro de los límites de precisión dada, los 
DÚMETOS Xy Y Xpt 


à n ? n 
Ejemplo. Hállese el valor de |/970 con un error inferior a 


0,001. Para k=3, la fórmula de aproximaciones (27) adquiere 
la forma 


a+2x 


E es) 


Apr = 


En nuestro caso a = 970. Hagamos x, = 10. De la fórmula (28) 
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se deduce que 
970+2-10% _ 2970 


a a T 
y n E299 _ 2910.60 _y 
ECO CS 


Vemos que dentro de los límites de la precisión dada los 
valores de x, y x, coinciden. Por ello, con un error inferior a 


0,001 tenemos 
V570 = 9,899. 


4 57 
MÉTODO DE ITERACIONES 


Todos los ejemplos expuestos más arriba son los casos par- 
ticulares de la aplicación de un método general para la resolu- 
ción de las ecuaciones. Este último se denomina método de ite- 
raciones o método de aproximaciones sucesivas. La esencia de este 
método consiste en lo siguiente. 

Una ecuación a ser resuelto, f(x) =0, se escribe en la forma 


x= Q(x). (29) 


Luego se elige la aproximación inicial X, y se sustituye en el 
segundo miembro de la ecuación (29). El valor obtenido de 
xa =p(x1) se toma por la segunda aproximación para la raiz. 
En general, si se conoce la aproximación Xm la aproximación 
posterior X,+, se halla según la fórmula 


Xati = P (Xa). 


Supongamos que después de realizar unas cuantas aproxima- 
ciones llegamos a la conclusión de que con el grado dado de 
precisión se cumple la igualdad x, = X„+1- Sabemos que Xn+1 = 
= (x). Esto significa que con el grado dado de precisión se 
cumple la igualdad x, = p(x,) es decir, que x, es un valor apro- 
ximado de la raiz de la ecuación x= (x). 

Por ejemplo, al resolver el problema sobre Aquiles y la tortuga, 
escribimos la ecuación 


10x = 1000 
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en la forma 
x 
x=10+ 10 
y buscåbamos las aproximaciones según la fórmula 
x, 
Xp+1 = 100 + w 


En el problema sobre la división con ayuda de una computa- 
dora escribimos la ecuación 


ax=b 
en la forma 
x=(1-a)x+b 
y buscábamos las aproximaciones por medio de la fórmula 
Xa+ =(1 ~ a) x, + b. 
Por fin, al extraer la raiz de k-ésima potencia, transformamos 


la ecuación 
x =a 
a la fòrmula 
ACES 

A 
después de lo cual buscábamos las aproximaciones según la 
fórmula 
a+(k- Dx 

kao 


Xati = 


Demos a conocer un ejemplo de la ecuación más compleja 
que se resuelve por el método de iteraciones. 


Ejemplo. Resuélvese la ecuación 
10x — 1 ~ cosx = 0 (30) 
con un error inferior a 0,001 


Escribamos la ecuación (30) en la forma 


1 + cosx 
x= i (31) 
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Elijamos alguna aproximación inicial, por ejemplo x; =0, y 
sustituyámosla en el segundo miembro de la igualdad (31) El valor 
obtenido 
L+c0s0 _ 

m 
lo tomaremos por la segunda aproximacion de la raiz buscada 
Sustituyendo el valor de x, en el segundo miembro de la igual- 
dad (31). obtendremos la tercera aproximación 


e 1450502 _ 14098 
z 0 10 


0.2 


Ha 


=0,198 


Luego hallamos 


= 145050198 
a= 10 = 0.198 
Vemos que con un error inferior a 0,001 se cumple la 1gual- 


1 + cosx 
dad xy = x,. Puesto que x¿= — io 2, esto significa que con 


un error inferior a 0,001 el número x,=0,198 es la raiz de la 
14c0sx 

10 

En relación con el método descrito de iteraciones surge toda 
una serie de preguntas. 

| ¿Siempre converge hacia un cierto número É la Sucesión 


Xp. Xp «+ Obtenida según el método de iteraciones? 
2. Si se verifica la igualdad lim x,=£%, ¿será el numero 3 la 
mx 


solución de la ecuacion x =p(x)? 

3» ¿Con cuál rapidez los números x;. .. x„ . Se aproximan 
a la raiz £ de la ecuación x = p(x)? 

Con mayor facilidad puede ser obtenida la respuesta a la segunda 
pregunta. Supongamos que los números Xy -s X «. SC apro- 
ximan al número £. Examinemos la igualdad x,.., = p(x). que nos 
proporciona la expresión para la aproximación posterior a través 
de la antecedente. Al aumentar n, el miembro primero se apro- 
xima a É, mientras que el segundo miembro se aproxima a 
p(E)* Por consiguiente, en limite obtendremos $ = tE) es decir, 
& es la raiz de la ecuación x = q(x) 


ecuación x 


+! Suponemos que «p(x) es una función continua 


3-71 
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La respuesta a la primera pregunta es negativa. Efectivamente, 
exploremos, por ejemplo, la ecuación 


x=10-2, 
Si hacemos x; = 1, obtendremos: 
x=8, x= 10 — 2... 


Con el aumento de n los números x, van creciendo y no 
tenden a ningún límite. Al mismo tiempo, si escribimos la ecua- 
ción en la forma x=Ig(x +2), el proceso de aproximaciones va 
a converger y, después de tres aproximaciones, obtenemos que 
x = 2,38, 

De este modo, la primera pregunta se la debe cambiar por 
otra que sigue: ¿cómo debe ser la función p(x) para que la 
sucesión de los números x1, Xz, ..., Xw .-. Converja? 

Antes de contestar a esta pregunta demos a conocer la inter- 
pretación geométrica del método de iteraciones. 


$8 ; 
SIGNIFICADO GEOMETRICO 
DEL METODO DE ITERACIONES 


Es evidente que la búsqueda de la raíz & de la ecuación 
x=0(x) significa nada más que la tentativa de la abscisa del 
punto M en el que se intersecan la curva y =«p(x) y la recta 
y=x. Sea x, (fig. 2) un cierto valor inicial. Entonces, el punto 
M, con coordenadas M, (X1, P(x1)) se ubica en la curva y =p (x). 
Tracemos por este punto una recta horizontal. Esta cortará la recta 
y=x en el punto N, (p(x1), p(x;)). Designemos y (x;) mediante 
xz. En este caso el punto N, tendrá las coordenadas N, (xz, x2). 
A continuación tracemos por el punto N, una recta vertical. 
Esta cortará la curva y=«Q(%) en el punto M2 (xa (x2)). 
Repitiendo este procedimiento, obtenemos el punto N; en la recta 
y =x con las coordenadas Nz(x3, x3), donde x, = Ọ (xz), y luego 
el punto M, en la curva y =(x) con las coordenadas Mj3(x3, 
(x3)), etc. Si el proceso de aproximaciones converge, los puntos 
My, Mz .... Mp ... van aproximándose al punto de inter- 
sección buscado. 

De esta manera, el significado geométrico del método de apro- 
ximaciones sucesivas consiste en lo que nos desplazamos al punto 
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buscado de la intersección de la curva con la recta a lo largo de 
una quebrada cuyos vértices se ubican alternativamente en la curva 
y en la recta, mientras que los lados tienen alternativamente las 
direcciones horizontal y vertical (fig. 2, a). 


Fig 2 


Si la curva y la recta están dispuestas tal como las muestra 
la fig. 2, a, la quebrada se parece a una escalera. En el caso 
representado por la fig. 2, b la línea quebrada es parecida a una 
espiral. 


Fig 3 


El proceso descrito de aproximaciones puede también divergir 
y no conducir a ningún resultado (igual que en el caso del problema 
sobre Aquiles y la tortuga). Gráficamente esto implica que los 
escalones de la escalera (o los eslabones de la espiral) se hacen 
cada vez mayores, por lo que los puntos My, ..-. Mm ... NO se 
aproximan al punto M sino se alejan de él (fig 3). 


ES 
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La diferencia entre las figs. 2 y 3 consiste en lo siguiente. 
Tracemos por el punto M, en el que se intersecan la recta y= x 
y la curva y= Q(x), una recta que forme con el eje de abscisas 
un ángulo de 135". Esta última recta junto con la recta y= x 
divide el plano en cuatro partes. Si la curva, que pasa en algún 
entorno del punto M, está dispuesta en los cuadrantes izquierdo 
y derecho del plano'y la aproximación inicial se toma en este 
entorno, entonces el proceso de iteraciones converge. Si en cambio, 
la curva está dispuesta en los cuadrantes superior e inferior del 
plano, el proceso de aproximaciones diverge. 

No obstante, para poder emplear esta regla, es preciso dibujar 
primero la gráfica de la función y = 6 (x), lo que no siempre parece 
conveniente, Por esta razón se debe establecer otro indicio de la 
convergencia del proceso de iteraciones que permitiría juzgar sobre 
la convergencia o la divergencia basándose en los cálculos anali- 
ticos y no de las construcciones geométricas. Esta condición será 
dada en el $ 10 Pero al principio demos a conocer el concepto de 
una aplicación contraida. 


$9 
APLICACIONES CONTRAIDAS 


Examinemos ta función y = (x), definida en el segmento [a, b]. 
A todo xo de este segmento le corresponde un punto yọ en el 
eje de ordenadas, a saber, el punto yo = p(xo). Para obtener este 
punto se debe trazar por el punto xy del eje de abscisas una 
recta vertical hasta que ésta corte la gráfica de la función 
y =Q(x), y luego trazar por el punto de intersección una recta 
horizontal hasta que se corten esta recta y el eje de ordenadas 
(fig. 4). De este modo, la función y=«Q(x) nos da la aplica- 
ción del segmento [a, b] en el eje de ordenadas. El conjunto 
de todos los puntos del eje de ordenadas correspondientes a los 
puntos del segmento [a, b], se denomina imagen de este segmento. 
Por ejemplo, la imagen del segmento [2,5] en la aplicación y = x? 
será el segmento [4,25], la imagen del segmento [ — 1,6] en la misma 
aplicación será el segmento [0,36] (sea al cargo del lector construir 
la gráfica de la función y= x°). Se puede demostrar que si la 
función y = (x) es continua en el segmento [a, b}, la imagen de 
este último será un segmento del eje de ordenadas. En el caso 
en que la función y =«p(x) sea, además, monótona creciente, la 
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imagén del segmento [a, b] se representará por el segmento 
Ri (DJ), y si sea monótona decreciente, por el segmento 
pb), ea] (ig. 5. 


ES 


Fig. 4 


En lugar de la aplicación del segmento [a, b] en el eje de 
ordenadas se puede considerar la aplicación suya en el eje de 


abscisas. Para ello, después de aplicarlo en el eje de ordenadas 
hagamos girar el eje mencionado 90” en el sentido horario. 
De resultas, los puntos del segmento [a, b) se aplicarán primero 
en ciertos puntos del eje de ordenadas y luego. en puntos del eje 
de abscisas. De este modo, la función p(x) determina la aplica- 
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ción del segmento [a, b] en el eje de abscisas. Esta aplicación 
la vamos a escribir así: x— ọ(x). Si p(x) es una función conti- 
nua, obtendremos, como resultado, algún segmento del eje de absci- 
sas. 

Puede ocurrir que la imagen [a,, b,] del segmento [a, b] 
sea una parte de [a, b]. Por ejemplo, en la aplicación y =/'x + 1 
el segmento [0,4] pasa a ser su propia parte [1,3]. Suele decirse 
en este caso que la aplicación «p(x) transforma el segmento [a, b) 
en la parte de éste. Si p(x) transforma el segmento [a, b] en la 
parte de él [a,, b,], entonces cualquier parte del segmento [a, b} 
se convierte en la parte del segmento [a,, b,]. En particular, 
el propio segmento [a,, b,] pasa, realizándose la aplicación p(x), 
en la parte de él mismo [a,, bz]. Del mismo modo, el segmento 
foa b,] se transforma en ła:misma aplicación a la parte propia 

az, by], etc. De resultas, obtendremos el sistema de segmentos 
[a, b], [as br), ... [ap b,], -.., cada uno de los cuales es una 
parte del segmento antecedente y de los segmentos tales que en la 
aplicación q (x) el segmento [a,+ 1, b,+1] sea la imágen del [a,, b,]. 


Por ejemplo, la aplicación x> 1 — z - 7 transforma el seg- 


mento [0,4] en ta parte de él ['/z, 5/6]. Aplicando esta apli- 
cación al segmento [*/2, */4] obtendremos el segmenta [?/s, **/17], 
etc. Cada segmento siguiente constituye una parte del antecedente. 

Se consideran dos casos eventuales: o bien existe el segmento 
[c, d], perteneciente a todos los segmentos [a,, b,], o bien estos 
segmentos tienen un solo punto común č. En el último caso se 
dice que el sistema de segmentos [an b,] se encoge al punto E. 

Más abajo enunciaremos la condición en la que se encoge 
el sistema de segmentos [a, b], [a,, by), ..., [ay b,], ... Para 
ello introduciremos un concepto importante de aplicación .con- 
traída, Una aplicación p(x) que transforma el segmento [a, b] en su 
parte [a,, bı], se llama contraída, si ella disminuye la distancia 
entre dos puntos cualesquiera de este segmento al menos en M 
veces, donde M > 1. Puesto que la distancia entre los puntos 
X1 Y xı mide |x2 —.X1l, la condición de contracción se puede 
enunciar así: 

Existe un número q tal que 0<gq<l y para cualesquiera dos 
puntos X,, x, del segmento [a, b] se cumple la desigualdad 


lo(x)— 0(x0| <q|x, —x,] (32) 
(aquí, q = 1/M). 
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La longitud de una parte arbitraria [c, d] del segmento [a, b] 
disminuye, en la aplicación contraída (x) por lo menos en 
M = 1/q veces. Efectivamente, sea [c,, d,] Ja imagen del segmento 
[c. d}. Entonces, cı y d, serán las imágenes de ciertos puntos 
x, y xz del segmento [c, d], 


c, = ọ(x;) di = p(x). 


Pero, en este caso 
ld, — cal = lp (x2) — 9l < q lx = xıl. 


Ya que los puntos x, y xz están dispuestos en el segmento 
c, d], la distancia entre ellos |xz — x, | es inferior a la longitud 
d —c| del segmento [c, d]. Quiere decir, 


|d -c | <alļd- c|. 


Nuestra afirmación queda demostrada. 

Ahora podemos enunciar la condición para que un sistema de 
segmentos [a,, bı), ..., [am ba], que se obtienen del segmento 
[a, b] aplicando sucesivamente la aplicación «p(x), se encoja. 

Si la aplicación y (x), que transforma el segmento [a, b] en una 
parte suya [a1, b], contrae dicho segmento, entonces el sistema de 
segmentos |a, bi), ..., [dn ba), ... se encgje a cierto punto E 
del segmento [a, b]. 

Efectivamente, siendo contraida la aplicación y (x), para todo n 
tenemos 


lbn — an| < glbri = an=: |. 

De modo igual 
|bn-1 — an- 1159 |ba-2 — a-21- 

Pero, en este caso 

lb, -a| < 4° [bn-2 ~ a-l- 
Repitiendo este razonamiento, obtenemos 

lba- a| <a |b- a]. 

Ya que 0 <q < Í, la sucesión de los números 4, 4°, .... q", -. 
tiende a cero, por lo cual las longitudes lo, - al de los segmenti 
[a,. b,] tienden a cero cuando n tiende a infinito. Pero, en esté 


caso no puede haber un segmento [c, d] que pertenezca a todos 
los segmentos [an b,]. Por consiguiente, el sistema de los segmentos 


[a, b}, [an bi] -.., [Qm ba), -- 
se encoje. 
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En conclusión examinemos las aplicaciones p (x) para las cuales 
la desigualdad (32) 


lo- p| <q) xh > 


donde 0 < q < 1, se cumple para cualquier par de números x; y Xz. 
Estas aplicaciones son contraídas en todo el eje numérico. Mos- 
tremos que en este caso existe un segmento que es contraido 
mediante la aplicación p(x). Dado que la condición (32) se cumple 
para cualesquiera dos puntos x, y X}, es suficiente demostrar que 
existe un segmento que es transformado en si mismo por medio 
de la aplicación p(x). Tomemos un número a cualquiera y hagámos 
b= (a). Elejamos el número q, < 1 tal que q < q1» 
Pongamos 


lb -a| 
1=4ı 


y mostreníos que el segmento [a — R, a + R} pasa, realizándose 
la aplicación ọ(x), en una parte de él mismo. En efecto, sea x 
un punto de este segmento. Entonces, x= a] <R. Debido a la 
desigualdad (32) encontramos que 


low- bl =|0() - pla] <q |x= al < aR. 
Pero, en este caso 
lp () — al = lp (x) —b+b=aj < 
<|o(x)-b]+]b-al <aR+|b=a|= 
>qR+(1—q)R=(1+9-q)R<R 


R 


Esto precisamente nos convence de que todo punto del segmento 
[a — R, a+ R] se transforma en la aplicación p(x) a un punto 
del mismo segmento, razón por la cual la aplicación p(x) contrae 
el segmento [a—R, a +R] > 


$ 10 ; 
APLICACIONES CONTRAÍDAS 
Y EL MÉTODO DE ITERACIONES 


Volvamos ahora al método de iteraciones. Este método se 
emplea para resolver las ecuaciones del tipo x = ọ (x) Si £ es la 
raíz de esta ecuación, £ = 9 (E), entonces, en la aplicación x p(x) 
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el punto E se queda en su lugar. De este modo, el problema 
de la resolución de la ecuación x= p(x) es equivalente al problema 
de la búsqueda de los puntos inmóviles de la aplicación p(x). 
Si la aplicación p(x) es contraida en el segmento [a, b], 
en este segmento siempre hay un punto inmóvil. Con el fin de 
convencernos de esto, tomemos una sucesión de los segmentos 


Ear biJ, [az ba), +, [tm ba), -+ 


que se obtienen del segmento [a, b) como resultado de efectuar 
sucesivamente la aplicación p(x). Ya que ọ(x) es la aplicación 
contraída del segmento [a, b], existe, por tanto, un punto único É, 
que es común para todos los segmentos [a,, b,]. Este es el punto 
inmóvil de la aplicación p(x) 

Efectivamente, la aplicación ọ(x) transforma todo segmento 
[aw b,] en su parte [a,,1. ba+:]. Por esto, la imagen «p(x) de 
todo punto c del segmento [dp ba] se ubica en su parte [an+ 1, 
Ba+1), y, con mayor razón, en el propio segmento [a,, b,]. Dado 
que el punto É pertenece a todos los segmentos [a,, b,], también 
su imagen ọ(Ẹ) debe pertenecer a todos estos segmentos. Pero, el 
único punto que pertenece a todos los segmentos [ap b,] es el 
punto É. Por consiguiente, p (E) =É, es decir, É es el punto in- 
móvil de la imagen p(x). 

Asi pues, para las aplicaciones que contraen el segmento [a, b], 
siempre existe un punto inmóvil ubicado en el mismo segmento, 
Es el punto único. En efecto, si hubiera otro punto inmóvil n, 
n = p(n), debería cumplirse la desigualdad 


In —8l=10(m) - e(8)] <qln—5l 


Como 0 < q < 1, esta desigualdad sólo puede tener lugar cuando 
|n- £| =0, es decir, cuando n = E. 

Podemos, ahora, enunciar la condición suficiente para la con- 
vergencia del proceso de iteraciones. 

Supongamos que la función (x) prefija la aplicación contraida 
del segmento [a, b]. Entonces, para cualquier punto xp de este 
segmento, la sucesión de múmeros Xo, Xi, Xz + s Xm o donde 
X„+ı = Ọ(x), converge hacia la raiz E de la ecuación x= Q(x), 
ubicadu en este segmento- 

Efectivamente, sea [ap b,}, ” = 1.2. ---, una sucesión. de los 
segmentos que se obtienen de [a, b] como resultado de efectuar 
sucesivamente la aplicación Q(x). Como el punto xo se ubica en 
el segmento [a, b], su imagen x, = P(x0) está cn el segmento 
[e bi], la imagen x, =p (x1) del punto x, se encuentra en el 
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segmento [az, bz], etc. Asi pues, para todo n, el punto x, se ubica 
en el segmento [am b,]. Puesto que las longitudes de los seg- 
mentos [am b,] tienden, con n creciente, a cero, la sucesión de 
los puntos Xp. .... Xp... tiende al punto común E de estos 
segmentos. 

El razonamiento realizado enseña que por el punto inicial xy 
puede tomarse cualquier punto del segmento [a, b). 

Aclaremos, ¿a qué velocidad tienden los puntos Xo, Xis <.. 


Xw -. al punto E? Como £= (E), para cualquier punto ¢ del 
segmento [a, b] tenemos 

lo()=El =19(0) -p| <qlc—85l. (33) 
Apliquemos ta desigualdad (33) a los puntos Xo, ..., Xp... Ya que 


Xa = P(x,- 1). entonces 
Lx —5l=1 0-1) El < q1X,-1 El. 
Pero, en este caso, para todo n tenemos 
la l< aln- El <q lxn-2 El <...<q"lxo —5l. 


De este modo, al crecer n, el error |x,—E| va disminuyendo al 
menos a la velocidad de una progresión geométrica con la razón q. 
Demos algunos ejemplos de la aplicación demostrada. 
Ejemplo 1. Puede aplicarse el método de iteraciones para 
resolver la ecuación 
1 


=372? (34) 


x 
En este caso 


=> 
PAE > 


Para cualesquiera x, y x, tenemos 
d epa 8 
4+ 44% 
lo al 
(DAD td 


Lo tx) — o (x1) 


be = xil- 


Pero, segùn la desigualdad entre el medio aritmético y el medio 
geometrico 
44+x? 


hl=> yx < Th 
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Por esta causa 
(4+ xå) + (443) 


Ixi + xl < lal +x sal 
2 2 2 2 2d 
aap E cay A y i STERKER 


Hemos demostrado que para cualesquiera x, y xz se cumple 
la desigualdad 
Xi +X 


UREA <i 


y, por tanto, 
1 
leixa) = pll <l- l 


Esto significa que la aplicación p(x) es contraída a lo largo de 
todo el eje. 

Ya sabemos que en este caso existe un segmento que en 
esta aplicación se transforma en si mismo. Para hallarlo, hagamos 
a =0, En la aplicación p (x) el punto a = 0 se transforma en el punto 

= 1/4. Luego, en nuestro caso q = '/g. Hagamos q, = 1/4 y desig- 
lb—al Eit 
=z. El =>». 3 

1% segmento 33 se 

transforma, realizándose la aplicación p(x), en sí mismo. Esto 


quiere decir que en este segmento existe un punto inmóvil, es decir, 
la raíz de la ecuación (34). Para hallar este punto, tomemos cual- 


nemos por R el número 


e £ T á 
quier punto del segmento [5 5) sea por ejemplo, xp =0. 


Empleando el método de iteraciones, obtenemos 
[i 
xı =7=025, 
x= l 
27025 N WE 
1 
A+ 024612 — T 
1 t 
“s = 44024637 — 40605 


qaar = 0,2461, 


= 0,2463, 


x= 


= 0,2463, 


Con un error inferior a 0,0001 tenemos xy = x4. Por consiguiente, 


la raíz de la ecuación (34), ubicada en el segmento [ 3] es igual 


a 0,2463 con un error inferior a 0,0001. Ya que la aplicación p(x) 
es contraída en todo el eje, la ecuación (34) no tiene otras raíces. 
Ejemplo 2. ¿Puede aplicarse el método de aproximaciones 
secesivas para resolver la ecuación 
3- 
x=l]+px 
en el segmento [— 1,8}? 


Ki 

Aquí, p(x)=1+|/x. Ya que p(—1)=0, (8) =3, entonces, 
la aplicación p(x) transforma el segmento [— 1,8] en si mismo. 
Sin embargo, la aplicación dada no es contraída en este segmento, 
puesto que, por ejemplo, para x, = 0,008, x, = 0,0008 


3 — 3, —— 
leixa) — 1x1) 1 = | V 0008 — )=0,008| = 0,4 > [x, = x, 1. 
Al demostar el hecho de que la aplicación en el ejemplo 1 fue 
contraida, empleamos la desigualdad /ab < “ 4 b Ahora deducirómos 
T 
0 A 


Fig. 6 


algunas desigualdades las cuales se usan con frecuencia para la demos- 
tración de que dicha aplicación es contraída. 
Demostremos que para x > 0 tiene lugar la desigualdad 


senx < x <tg(x) (35) 
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Notemos, para ello, que en la fig. 6 el área Soag del sector OAB con el 
ángulo central x está comprendida entre tas áreas de los triángulos 
OAB y OAT: 


S 2048 < Siecror 04B < S 2047 
Pero, 


R? sen x Rè tgx 
Saos = > Saorr= z 


(R es el radio del circulo). Mientras tanto, el área de) sector OAB es 
igual a en (medimos el ángulo en radianes). Por esto, 
Riseax Rx Ritgx 
2 2 2 
Al reducir por a llegamos a la desigualdad (35). De la desigualdad (35) 


se desprende que para O < x < 1, tenemos 


x < arcsenx, 


y para x > 0, tenemos 
Xx > arctg x. 


Indiquemos también las desigualdades 
é>i+x x>0 
mä +x) <x, 0<x<1l 


cuya demostración es un poco más complicada. 


Ejemplo 3, Aclárese, si puede resolverse mediante el método 
de iteraciones la ecuación 


x=1+ Jaragx (36) 
Dado que para todo valor de x se tiene 1 +3 arctg x > 0, entonces 
la ecuación sólo puede tener las raices positivas. Tenemos: 


o()=1+ jarag xX. 7 


Por esto, 


lo(x) — e(x1)l= |( + Jara) = ( + Jara +) 


larcig xz — arctg xı |. 


ES 
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Pero, para x,>0, x2>0 


X2 = X; 
arctg xz — arcig x; mate ier 
X1X2 
y, por lo tanto, 
1 X= 1| xx 1 
- => |arctg ——_L| <A Rs xl. 
lo (x2) = Q (x1) = 5| arctg Er Dl rr a L 


Esto significa que la aplicación (37) es contraída en el semieje [0, 00]. 
Esta aplicación transforma el segmento [o, va] en su parte 
[1, 1 + 1/6]. Por esta razón, existe una raíz única de la ecuación (36), 
ubicada en el segmento [1, 1 + 1/6]. Para hallar esta raíz, haga- 
mos x, = 1. En este caso, 

n 


3 = 1,39, 


x= + zorctgl =1+ 
x=1 + jarag 1,39 = 1,474, 


Gas garag 1,474 = 1,487, 


xs = 1 + arctg 1,487 = 1,489, 
sje zaag 1,489 = 1,490, 


x=l+ jarig 1,490 = 1,490. 


Vemos que la igualdad xe = x, = 1,490 se cumple con un error 
inferior a 0,001. Quiere decir que con esta precisión la raíz de 
nuestra ecuación es igual a 1,490. Ya que la aplicación p(x) es 
contraída en todo el semieje 0 < x< œ, resulta que la ecuación 
(36) no tiene otras raíces. 

Ocurre con frecuencia que la ecuación x =p(x), la cual no se 
resuelve por el método de iteraciones, puede transformarse a una 
ecuación del tipo que admite este método. Tomemos, por ejemplo, 
la ecuación 

x=x -2 (38) 
Puesto que aqui 


ọ(l)= -1 <1, 9(2)=6>2, 
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esta ecuación tiene una raiz en el segmento [!, 2]. Pero, la apli- 
cación x? — 2 no es contraída en este segmento, ya que no aplica 
este segmento en su parte. Escribamos la ecuación (38) en la forma 


3 
x=hx+2. 
Aquí, Y (x) = ix +2 y, por tanto, 
A US 
Ie) = vidls ||% +2- px +2|= 


XX 


7 F) 7 x 
| + Y + Ax + 2) +D hix +2) 


En el segmento [1, 2) tenemos x, > 1. x,>1. Por consiguiente, 
en este segmento 


Iyad- yels 


315 
Hemos demostrado, pues, que la aplicación y(x) es contraida en 
el segmento [1, 2]. Hagamos x; = 1 y apliquemos el método de 


iteraciones. 
Obtendremos 


lx xx) 


x=13=1442, 

xa = 13,442 = 1,510, 
xa = [3310 = 1,520, 
xs = | 3520 = 1,521, 
xe = 3521 = 1,521. 


Así pues, la raiz de la ecuación (38), ubicada en el segmento 
[1, 2] es, con un error inferior a 0,0001, el número 1,521, Esta 
ecuación no tiene otras raices. 

Vemos que la transformación acertada de la ecuación de partida 
la llevó a la de otro tipo que admite usar el método de iteraciones. 

El criterio de convergencia del método de iteraciones, descrito 
aquí, no es muy cómodo para el empleo, dado que requiere la 
demostración de las desigualdades bastante complejas. A conti- 
nuación {en el $ 21) vamos a conocer un corolario de este criterio, ` 
que simplifica considerablemente el procedimiento de demostrar la 
convergencia del proceso de iteraciones 
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h $11 
METODO DE CUERDAS 


El método de iteraciones es uno de los métodos más generales 
de la'resolución aproximada de ecuaciones. Muchos otros métodos 
para la resolución aproximada de ecuaciones representan casos 
particulares del método de iteraciones. Ahora, demos a conocer 
uno de estos métodos, llamado método de cuerdas. 

Supongamos que se necesita resolver la ecuación Ja =0 
Este problema es equivalente al de búsqueda de los puntos, en los 
cuales la curva de la función y= f(x) corta el eje de abscisas. 


Fig. 7 


Supongamos también que la función f(x) es continua y en los puntos 
a y b tiene valores de signos contrarios. Entonces, al menos en un 
punto del segmento [a, b] esta función se anula. En otras pala- 
bras, al menos en un punto del segmento [a, b] el gráfico de la 
función y = f(x} corta el eje de abscisas. En general, pueden haber 
unos cuantos de tales puntos (fig. 7). No obstante, si la función 
y = f(x) es monótona en el segmento [a, b] y tiene en sus extremos 
los valores de signos opuestos, el gráfico de la función en este 
segmento se interseca con el eje de abscisas sólo en un único 
punto č. Para hallar aproximadamente este punto, sustituyamos 
en el segmento [a, b] el arco de la curva y=f(x) por la cuerda 
correspondiente MN y encontraremos el punto de intersección de 
esta cuerda con el eje Ox (fig. 8). 

Examinemos con este fin los triángulos semejantes MMT y 
NN,T. De la semejanza de estos triángulos se deduce que 
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Y SL Pero, de la fig. 8 se ve que M¡T=a,-—a, 


TN, =b-a,, MM, = — fla) y NiN = f(b) donde mediante a, 
está designada la abscisa del punto de intersección de la cuerda 


Y N 


Fig. 8 


MN con el eje Ox. Por esta causa 
aa b-a, 
-A TÖ 
Al resolver esta ecuación, encontramos que 
am ab) — bfla) ] 
fib) — fa) 


Esta igualdad puede escribirse también en la forma 
b-a 
e 1075=1 69 
*o bien k 
b-a 
00 ee 


(compruébese, reduciendo los segundos miembros en las fórmulas 
(39) y (40) a un mismo denominador). 

El número a, es precisamente el valor aproximado de la raíz 
de la ecuación f(x) =0, que se encuentra entre los puntos a y b. 


4-71 
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Como los signos de los números fía) y f(b} son opuestos, 
entonces, o bien el signo de fía) o bien el de f(b) se diferencian 
del signo de f(a,). Si en los puntos a y a, la función f(x) tiene 
signos contrarios, se aplica al segmento [a, a,] la fórmula (39) 
y se obtiene la aproximación apan para la raíz buscada: 


far) 


ON en 


Si la función f(x) asume los valores de signos opuestos en los 
puntos a, y b, se aplica la fórmula (40) al segmento [a,, b] y se 
supone 


a =m -Sd ET (42) 


Al hallar el valor az, se aplica la fórmula (39) al segmento 
[a, a2] (o, correspondientemente, la fórmula (40) al segmento [a,, b]) 
y se halla la siguiente aproximación az. En general, una vez hallada 
ba aproximación a,, se busca la siguiente aproximación según la 
fórmula 


a — a 

n41 = 4, — Sa.) Ta) (43) 
o, según la fórmula 

ansi = n= Ma Ha)" (44) 


Hemos obtenido dos fórmulas (43) y (44). Veamos, ahora, en 
qué caso debe ser usada una y otra fórmula. Supongamos que la 
curva gira su concavidad hacia la parte positiva del eje de ordenadas. 
En este caso se deben unir los puntos de la curva con aquel 
extremo M O N, en el que la función es positiva. Si, en cambio, la 
curva gira su concavidad hacia la parte negativa del eje de orde- 
nadas, los puntos de la curva se unen con aquel extremo, en el 
cual la función es negativa. En la figura 9 aparecen los diferentes 
casos posibles. Estos dibujos hacen evidente geométricamente la afir- 
mación siguiente. 

Supongamos que en el segmento [a, b] la función f(x) es continua, 
monótona, gira su concavidad hacía la dirección constante y admite 
en los extremos del segmento los valores de signos opuestos. En este 
caso, con la elección adecuada de la fórmula de aproximaciones el 
método de cuerdas da una sucesión de puntos convergente hacia la 
raíz de la ecuación f(x) = 


31 


En el caso de que la elección de la fórmula no sea acertada, 
el método de cuerdas puede llevar a que el punto az resulte 
fuera de los límites del segmento [a, b). Esto se ilustra en la fig. 10. 

El método descrito de cuerdas es un caso particular del método 


o) d) 


Fig. 9 


—_ == 


de iteraciones. Supongamos que la función f(x) no se reduce a cero, 
cuando x =a. Entonces, la ecuación f(x)=0 es equivalente a la 
ecuación 


x-a 
SOS Gi 
En efecto, si f(£) = 0, entonces 
ÉS Ea 
E= -SOENE (46) 


Inversamente, si >a y se cumple la igualdad (46), entonces 


JE) =0. 


“e 
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Pero, la ecuación (45) tiene por expresión x= ọ(x) donde 


x-a af(x) — xfa) 
=x- q P a e r e 
CIAO TO 
Hagamos xp =b y apliquemos el método de iteraciones. Obten- 
dremos la misma sucesión de números Bis Az, ooy Amy ossis QUE 
tuvimos al emplear el método de cuerdas: 


Fig. 10 
Sati g e e a a 


Como ejemplo, resolvamos, por el método de cuerdas, la ecuación 
34+3x-1=0, (47) 


Aquí, f(x) =x? + 3x— 1. Dado que f(0) = —1, f(1)=3, entonces 
la ecuación (47) tiene por lo menos una raíz en el segmento [0, 1. 
Si dibujamos el gráfico de la función y =x? +3x-— 1, veremos 
que en el segmento [0, 1] él gira su concavidad hacia la parte 
positiva del eje de ordenadas. Por esta razón aplicamos la fórmula 
(39). Según la fórmula (39), la primera aproximación para esta 
raíz es el número 


b~a 
bi 
La segunda aproximación se halla por medio de la fórmula 
b- xı 1—0, d 
=b- fib- = 1-3 Mm 
xXx =b-fi e = FA 3 37023 0,31 


A continuación, 


1-031 _ 
x= 1-3 6% = 031 
o 10319 _ 
13 oro = 032, 
1 -0,322 
xs=1- 3y 000057 0322 


De este modo, con un error inferior a 0,001, la raiz de la ecua- 
ción, ubicada en el segmento 0,1, es igual a 0,322. 


z yu 
MÉTODO DE CUERDAS 
PERFECCIONADO 


Si el método de cuerdas converge, la velocidad de la conver- 
gencia será la misma que es propia para el método de iteraciones: 
el error en el valor de la raiz disminuye como una progresión 
geométrica. El perfeccionamiento del método de cuerdas da la 
convergencia mucho más rápida. En el método de cuerdas habitual 
se usa, en cada paso, uno de los extremos del segmento [a, b] 
y la última aproximación obtenida. En lugar de esto pueden ser 
empleadas dos últimas aproximaciones, pues, ellos son, comúnmente, 
más cercanas a la raíz buscada que los extremos del segmento 
[a, b}. 

La fórmula, en la cual se utilizan dos últimas aproximaciones, 
tiene la forma siguiente (fig. 11, a): 


a n= Ona 

Or IO pa fa "a 
Con ello, a, se calcula por medio de la fórmula (39) y az, 
mediante las fórmulas (41) ó (42), según sean los signos de f(a), 
fib), f(a,): si fla) < 0, f(b) > 0, entonces para f(a) < 0 se toma la 

fórmula (42), y cuando f(a1) > 0, se elige la fórmula (41). 
Si, por casualidad, llegamos a que el punto ay, calculado por 
medio de la fórmula (48), está fuera de los limites del segmento 
[a, b], entonces, realizando el paso siguiente, en lugar de este punto 
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se debe tomar el extremo del segmento [a, b] que sea más próximo 
al punto citado (fig. 11, b). 

Resulta pues, que la convergencia del método perfeccionado de 
cuerdas es más rápida que la del método ordinario. A' saber, 


a) 6) 
Fig. 11 


si É es una raíz de la ecuación f(x) =0, entonces 


lapsi El < Clan — El", (49) 


LA 1,618. 


A título de ejemplo usaremos este método para resolver la 
misma ecuación 


donde 


x +3x-1=0, 


que antes resolviamos con ayuda del método de cuerdas. Las 
primeras aproximaciones, a, = 0,25 y a, = 0,31, son las mismas que 
en el método habitual de cuerdas. 

La siguiente aproximación calculamos por medio de la fórmula 


0,31 —0,25 
2040702347 0,3223. 


a— a, 
=a- Foi Aa =0,31+0,040 
dd e oE 
Obtenemos que /(0,3223) = 0,0004, Está claro que x= 0,3223 
nos da la raiz buscada con un error inferior a 0,0001. 
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$13 
DERIVADA DE UN POLINOMIO 


Al resolver la ecuación f(x) =0 por el método de iteraciones, 
mucho depende del método por el cual la ecuación citada fue redu- 
cida a la forma x= ọ(x) En muchas circunstancias resulta más 
conveniente el método propuesto por Newton. Este método se basa 
en el concepto de la derivada. En este párrafo vamos a relatar 
qué es la derivada de un polinomio. Esto nos permitirá aplicar 
el método de Newton a la resolución de Jas ecuaciones algebraicas, 
es decir de las ecuaciones del tipo 

aox + aX" +... +a =O. (50) 
Sea 
f)=ax*+ax H+. +4 


— cierto polinomio. Examinemos el polinomio f(x + a), es decir, la 
expresión E 

alta + alx +a +.. H yo (51) 
Si abrimos paréntesis en la expresión (51), en algunos sumandos « 
faltará, en otros sumandos figurará en el primer grado, en los 
terceros en el segundo grado, etc. Reunamos los términos que 


contienen œ de un mismo grado, Entonces, el polinomio f(x + a) 
tendrá por expresión 


fra =f tA daras hna (59 


(ya que el grado del polinomio f(x) es igual a k, el grado supe- 
rior de a en el desarrollo (52) también es igual a k). Es evidente, 
además, que fo(x), ..., f(x) también son polinomios de x. 


Ejemplo. Sea 
f(x) =2x?— 3x? 4 6x— 1. 
Entonces, 
fix +a) =2Ux +a}? — Fx +a)? + 6(x+ 1 = 
= 2 + 3xta + 3xa? +a?) — 3 + 2x0 + a?) + 6(x +a) 1= 
= (2x? — 3x? + 6x — 1) + (6x? — 6x + 6Ja + (6x — Ya? + 200. 
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Por lo tanto, en este caso 

fo) =2x? — 3x? + 6x — 1, 

$0) = 6x? — 6x +6, 

hx) =6x — 3, 

Six) =2. 
Vemos que el sumando f(x) coincide con f(x). Esta coincidencia 
no es casual. Si en la igualdad (52) ponemos a =0, obtendremos 
que fix) = falx). . 

Exploremos, ahora, el siguiente sumando f, (x)a. El coeficiente de 
a, esto es, el polinomio f,(x), se denomina derivada del polinomio 
fix). Por ejemplo, la derivada del polinomio 2x? — 3x? + 6x — 1 
tiene por expresión 6x? — 6x — 6, La derivada del polinomio f(x) 
se designa mediante el símbolo f(x). 

Así pues, la derivada f' (x) del polinomio f(x) se llama el coefi- 
ciente de œ en el desarrollo del polinomio f(x + a) segisn los grados 
de a. 

Al emplear las designaciones introducidas, la fórmula (52) puede 
escribirse en la forma 

fx+o)=flx)+/f (a+... (53) 

Los puntos en esta fórmula denotan los términos que contienen 

aê, 00%, ..., al. Por ejemplo, 
2(x + a)? — 3(x +a)? + 6(x + a) = 1 = 2x? — 3x? + 6x1 + 
+ (6x? — 6x + Oa +... 

Hemos introducido el concepto de. la derivada de un poli- 
nomio f(x) Señalemos, ahora, cómo se calcula esta derivada. 
Examinemos con este fin el polinomio 


Jx+0)=a (x+ aF + a(x + a'H... + ao (X H 0) + oe 


Sustituyendo cada sumando según la fórmula (x + a" = x" + 
+ mx" 'a +... (véase $ 6), obtendremos 


fix +a) = abt tkt a+...) 
+0 [71 + (k— Dada] onn H aa l H a) H a= 
=a +a 0 7 ++ 
+a [ka + (k— I)a? +... +A]... 
Comparemos esta igualdad con la (53): 
fr = M0 +a +... 
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Obtenemos la siguiente afirmación: 
La derivada del polinomio 
arta +... + 1 A a (54) 
tiene por expresión 
f'(x) =kag 7! + (k — I)a? +... + 0-1 (55) 
Por ejemplo, la derivada del polinomio ý 
f(x) = 6x? + 8x? -3x* —1 


2 


f(x) = 42x* + 24x? — 6x. 


, 514 
MÉTODO DE NEWTON PARA 
LA RESOLUCIÓN APROXIMADA 
DE LAS ECUACIONES ALGEBRAICAS 


Volvamos a la resolución aproximada de las ecuaciones alge- 
braicas, Sea dada la ecuación 


aox +a? +... + a =O. (56) 


Supongamos que de una u otra manera fue hallado el valor 
aproximado x, de la raíz de esta ecuación y señalemos cómo 
se puede hallar el valor de esta raíz más exacto. Designemos 
con a, el error del valor xy, es decir, supongamos que x, + a, es la 
raíz de la ecuación (56). En este caso tiene lugar la igualdad 


ax +0) tala +) t+...+09=0. (57) 
En otras palabras, 
Sx +01) =0, 
donde con f(x) está designado el polinomio 
E t., H ay 
Pero, según la fórmula (31), tenemos 
Si +a) ES) AA Ae. 


donde por medio de los puntos están designados los términos 
que contienen af, ..., af. De este modo, para determinar %, 
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nos disponemos de la ecuación 
fx +a) = f(x) + aS (X) +.. (58) 


Si la aproximación inicial x, fue suficientemente acertada, el error 
a, de esta aproximación es pequeño. En este caso, los términos 
en la igualdad (58), designados con puntos, serán pequeños en com- 
paración con a,. Despreciando estos términos, obtendremos la 
igualdad aproximada para determinar a, 


Fx) + af (x) 0. (59) 

De esta igualdad se deduce que 
(60) 
Por eso, el valor mejorado x, para la raíz de nuestra ecuación 
se define por medio de la fórmula 
SA 

. 61 
SA) e 
Después de esto, puede ser mejorada de nuevo la aproximación 


hallada. La tercera aproximación para la raíz se da por medio 
de la fórmula 
JA) 


Y 

En general, una vez determinada la n-ésima aproximación x, de 
la raíz buscada, la siguiente aproximación se da por medio de 
la fórmula 


aA=x- 


Xy" Xz 


_ IG 
S (xa) 
En la forma desarrollada esta fórmula se escribe así: 
aa +... + 1% + ay 


ari 


Si, en los límites de la precisión prefijada, los valores aproxi- 
mados X, y X„+ı Coinciden, nuestro proceso (en los limites de la 
precisión prefijada) se considera terminado, y el valor buscado 
de la raíz está determinado. 

El método descrito de la resolución de las ecuaciones se debe a 
Newton, famoso matemático inglés. 

El método de Newton está estrechamente ligado con el método de 
iteraciones. A saber, si las funciones y = f(x) e y = f'(x) no tienen 


Xati = Xp 


(62) 
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raices comunes, la ecuación f(x) =0 es equivalente a la ecuación 
xxi, (64) 


Al aplicar a esta ecuación el método de iteraciones, obtendremos 
una sucesión de números Xy, Xz, ..., Xm ligados con la misma 
correlación 

ES) 


TE" 

que teníamos, al operar con el método de Newton. En otras 

palabras el método de Newton consiste en lo que la ecuación 

f(x) =0 se escribe en la forma (64) y a él se le aplica el método 

de iteraciones. 

Ejemplo. Resuélvase por el método de Newton, la ecuación 
*-3x-5=0 


Xari Xp 


(65) 


con un error inferior a 0,001, tomando por la primera aproxi- 
mación xı =3. 
Como la derivada del polinomio 


fa) =x -3x-5 
es el polinomio 


S0)=3-3, 
entonces, la fórmula (62) tomará, en este caso, la forma 
EITE AA 
AA 
Por esta causa, 
27-9-5 13 
x=3- RF “332 
14,89 — 7,38 — 5 e 
Xa = 2,46 — gea A 2,46 — 0,165 = 2,295 
12,088 — 6,885 — 5 
M4 = 2295 - — isg 3 "2095 0.016 = 2,279, 
E 11,837 — 6,807 -5 _ 
200-221 


Vemos que la igualdad 
X4 =Xs 
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se cumple con un error inferior a 0,001. Por esta causa, la raiz 
de la ecuación x?—3x-—>5=0 es igual a 2,279 con un error 
inferior a 0,001. 

El método del cálculo aproximado de raices, expuesto en el 
$ 6, es un caso particular del método de Newton. Efectivamente, 


la búsqueda de Va no es otra cosa que la resolución de la ecua- 
ción 
x*x—-a=0, 


Pero, la derivada del polinomio x*—a es igual a kx*”?, y, 
por tanto, para la ecuación 


*=a=0 
la fórmula (62) tendrá por expresión 


zd B- a a+ík-i)xŻ 
Xati = Xp ys DMA» e 


Esta gs la fórmula, según la cual se calculan las aproximaciones 
para \ a. 

Indiquemos la siguiente diferencia esencial que existe entre las 
soluciones de la ecuación x*—a=0 y de la ecuación algebraica 
general 


axt +a č! +... + a =O. 


Para la ecuación x* — a = 0, la elección de la aproximación inicial 
x, no tenía mucha importancia. Sea cual fuera el valor de xy, 
gl_haber realizado unos cuantos pasos, llegaremos al valor de 
pa con la precisión prefijada. No es así, cuando resolvemos la 
ecuación (56). En este caso, algunos valores iniciales llevan a ciertas 
raíces de la ecuación, otros valores iniciales conducen a las raíces 
diferentes, y ciertos valores iniciales no llevan a ningún valor 
determinado de la raíz: la sucesión de números Xi, Xz, +-+, Xp +... 
calculada por medio de la fórmula (62), no tiende a ningún 
límite determinado, es decir, es divergente. 
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$ 15 s 
SIGNIFICADO GEOMÉTRICO 
DE LA DERIVADA 


El método de Newton lo expusimos sólo para las-ecuaciones 
algebraicas. Con el fin de extenderlo a las ecuaciones de forma 
arbitraria, es necesario generalizar el concepto de la derivada, 
introduciéndolo para cualesquiera funciones. Con este objeto 


aclaremos el significado geométrico de la derivada. Examinemos 
el gráfico del polinomio 


y= agt + axt) +... a 


y escojamos en este gráfico dos puntos M y N (fig. 12). 
Supongamos que la abscisa del punto M es x, y en el punto 
N la abscisa es igual a x + a. Entonces, las ordenadas en los puntos 
M y N se definen, respectivamente, por las expresiones 


SO) = apt + axt n 


fix +a) = alx +a + a(x tal H.. H de 
Tracemos por los puntos M y N. una secante y calculemos el 
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Coeficiente angular k,.¿*? de esta secante. Del dibujo se ve que 
TN 
tey =—. 
ev= T 


Pero, el segmento M Tes igual a la diferencia entre las abscisas de los 
puntos M y N, y, por tanto, 


MT=(x+0)-—x=0, 
El segmento TN es igual a la diferencia entre las ordenadas de estos 
puntos y, por tanto, 
TN f(x + a) — fix). 
Por consiguiente, 
TN E — fix] 
y= ipa 14 sœ 
Pero, según la fórmula (53), 
Fix + =f YURI 


donde con los puntos están designados los términos que contienen 
C 
Por esta causa 


tgy = aft.. 


afix) +... 
a 
donde con los puntos están designados esta vez los sumandos 
que contienen a, x?, ... 
Asi pues, el coeficiente angular de la secante MN se expresa por 
medio de la fórmula 


Ku =t8Y=J (+... (66) 
Vamos ahora a disminuir el valor de a. La secante MN girará, 
en este caso, alrededor del punto M. En límite, para x=0, la 


secante se transformará en una tangente a la curva en el punto M. 
La fig. 13 ilustra las posiciones de la secante para «=l; 
E, f 


a 


* Se llama coeficiente angular de una recta a la tan- 
gente del ángulo que forma esta recta con el sentido positivo del eje Ox. 
Por ejemplo, si la recta forma con_el eje Ox el ángulo de 60”, el coeficiente 


angular de la recta es igual a |/3. 
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Pero, todos los términos, designados en la fórmula (66) con 
puntos, se anulan cuando a =0. Por eso, el coeficiente angular 
de la tangente al gráfico del polinomio y= fix) en el punto de 
abscisa x se expresa por la fórmula 


tee ES). (67) 


Así pues, la derivada del polinomio fix) es igual al coeficiente 
angular de la tangente al gráfico de este polinomio, trazado en el 
punto de abscisa x. 

Ejemplo. Hállese el ángulo que forma con el eje Ox la recta 
tangente al gráfico del polinomio 


f= -4x+51+4 
trazado en el punto de abscisa x = 2. 
Dado que f'(x) = 3x? — 8x +5, entonces f (2) = 1. Por consi- 


guiente, tg y = 1, de lo que se deduce que el ángulo de inclinación p 
es igual a 45°. n 
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$ 16 i 
SIGNIFICADO GEOMÉTRICO 
DEL MÉTODO DE NEWTON 


Ahọra podemos aclarar el significado geométrico del método 
de Newton para la resolución aproximada de las ecuaciones alge- 
braicas. Supongamos que se necesita resolver la ecuación f(x) = 0, 
donde f(x) es un cierto polinomio. Desde el punto de vista geo- 
métrico es un problema en que se buscan los puntos de intersección 
del gráfico de la función y=f(x) con el eje Ox, es decir, los 
puntos donde y =Q. 


y 


Fig. 14 


Supongamos que el valor aproximado x, de la raíz de esta 
ecuación ya está conocido. Tracemos por el punto N de abscisa 
xı una tangente a la curva y= f(x). Siendo acertada la elección 
del valor x,, el punto T de la intersección de la tangente con el 
eje Ox será más próximo al punto en que se cortan la curva 
y=f(x) y el eje Ox que el punto M (fig. 14). 

Para hallar la abscisa x, del punto T, examinemos el triángulo 
TMN. El cateto MN de este triángulo no es otra cosa que el valor 
de la función y=f(x) en el punto x,, es decir, MN = f(x,). 
Mientras tanto, el cateto TM es igual a xı — x De aquí se 
deduce que la tangente del ángulo q, que forma la recta tangente 
con el eje Ox se expresa por medio de la fórmula 


Se) 


X= Xa 


1891 = (68) 
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De la igualdad (68) se desprende que 
+. LEN) 
X2=X1 Bs (69) 
Pero tg, es el coeficiente angular de la tangente a la curva 
y= f(x), trazada por el punto de abscisa x,. Por esta razón, en 
virtud del significado geométrico de la derivada, tg 9, = f (x1). 
Esto significa que la fórmula (69) puede ser escrita así: 
SA) 
SUxp) 
Por consiguiente, hemos hallado la segunda aproximación para 
la raíz buscada. Tracemos, ahora, una tangente a la curva y = f(x) 


a=Xau- 


y 


Fig 15 


por el punto, donde la abscisa es xı. La abscisa del punto de 
intersección de esta tangente con el eje Ox se da por medio de la 
formula 

fx) 


O 

En general, una vez hallada la aproximación x,, entonces, 
para obtener la aproximación siguiente x„+ı se debe trazar una 
tangente a la curva y =f(x) en el punto de abscisa x,. La abscisa 
del punto de intersección de esta tangente con el eje Ox nos da X,+ 1- 


*x3=X 


5-71 
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La formula para calcular x,., €s: 


_ 10) 


(70) 


Xara = Xp 


S nd 


o, lo que es lo mismo, 
(10) 


donde q, es el ángulo de inclinación de la tangente a la curva 
y= fix) en el punto de abscisa x, Esta fórmula coincide con la 
fórmula (62) del método de Newton. Pusimos en claro, pues, 
el significado geométrico del método de Newton. El consiste en que 
el arco de la curva y= f(x) se sustituye por la tangente a esta 
curva. Por eso, el método de Newton se llama también método de las 
tangentes. En la fig. 15 se ve como los puntos Xy, Xz -+s Xm o 
que se obtienen por el método de Newton, van aproximándose 
al punto £, donde se cortan la curva y =f(x) y el eje Ox. 


$1 
DERIVADAS DE LAS FUNCIONES 
CUALESQUIERA 


El significado geométrico expuesto del método de Newton permite 
generalizarlo sin dificultad alguna para cualesquiera ecuaciones 
del tipo f(x)=0, donde f(x) ya puede no ser un polinomio. 
A saber, para hallar la solución de esta ecuación, tomemos algún 
valor aproximado x, de la raíz. Tracemos una tangente a la curva 
y = f(x) por el punto de abscisa x, y designemos con x, el punto 
donde la tangente corta el eje Ox. Por el punto de abscisa x, 
tracemos de nuevo una tangente a la curva y = f(x), etc. Al igual 
que en el caso en que f(x) sea un polinomio, es fácil establecer 
que 
_ (5) 

TN 


donde tg, es el coeficiente angular de la tangente a la curva 
y =/f(x) en el punto de abscisa Xy. 

Por el momento ta fórmula (71) no sirve para los cálculos, 
puesto que todavía no sabemos determinar tg p,. Por consiguiente. 
debemos aprender a calcular los coeficientes angulares de las tan- 


Xari Na 


(71) 
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gentes a los gráficos de cualesquiera funciones y = f(x) (y no sólo a los 
gráficos de los polinomios). Hallemos primero el coeficiente angular 
de una secante. Sea M un punto en el gráfico de la función 
y=f(x), y MN, una secante que pasa por este punto, Razonando 
igual que en el caso de los polinomios, obtenemos que el coefi- 
ciente angular de la secante se expresa mediante la fórmula 


Íx+ af) 


ke = 18 y = a 


(12) 
donde x es la abscisa del punto M, y x + a, la abscisa del punto N. 
Si hacemos disminuir a, la secante va a girar alrededor del punto 
M y tender a ocupar la posición de la tangente a la curva 
y = f(x) en este punto (véase fig. 12). Por esto, podemos escribir que 


Six +a) - fix) 


ka= 189 = lim EE 


(73) 
El limite en el segundo miembro de esta igualdad lo llamaremos 
derivada de la función f(x) y lo designaremos por el símbolo f fx), es 
decir, hagamos 


fí)= limeto A, aa 


Podemos, ahora, escribir la igualdad (73) en la forma 
ku = tgo = f'(x). (75) 


De este modo, también para las funciones cualesquiera (y no sólo 
para los polinomios), la derivada de la función f(x) en cierto punto 
es igual al coeficiente angular de la tangente a la curva y= f(x) 
en este punto”. 
Puesto que tg 9, =f’ (x,), la fórmula (71) puede ser escrita en 
la forma 
F(x) 


Ld 


Esta fórmula coincide con la fórmula (62). De este modo, el 
método de Newton queda generalizado para todas las ecuaciones del 
tipo f(x) = 0. 


Xara = Xan — (76) 


* Si por algun punto de abscisa x no puede ser tra- 
zada una tangente al gráfico de Ja función + = fiv) (por ejemplo, s: en este 
punto el gráfico sufre wn viraje), entonces la funcion y= fix) no tiene 
derivada en dicho punto 
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518 
CALCULO 
DE LAS DERIVADAS 


Hemos visto en el párrafo anterior que para hallar el coefi- 
ciente angular de la tangente a la curva y = f(x), se debe calcular el 
límite 

J'~) = lim 


Li+ a) -f 
2-0 a ` 

Este cálculo, en general, es bastante dificultoso. Pero, este 
límite está calculado para muchos casos importantes. En otras 
palabras, para las funciones que se encuentran con mayor frecuencia, 
sus derivadas están conocidas. Abajo vienen dadas derivadas más 
usadas: 


liak OEE. 
1. (ay = 0. 7. (ctgax) T 
2. (1) Aa. 
3. (a) = alna. 3. (log,x) = xlna ` 
4. (senax) = acosax. 9. (arcsenax) = K : 
S. (cosax) = — asenax. Yi -a 

j 10. (arctgaxY = — rr: 

6. (tgax) = TES gax) Teea 


(En las fórmulas 3 y 8 mediante In a está designado el logaritmo 
que tiene por base el número e = 2,71828..., o sea, el logaritmo 
natural). Señalemos que en la fórmula 2 el exponente k puede ser 
no sólo un número natural, sino también cualquier número real. Por 
ejemplo, 


Las fórmulas t— 10 son aún insuficientes para el cálculo de las 
derivadas de unas funciones cualesquiera. Sin embargo, si la función 
f(x) se compone, mediante las operaciones aritméticas, de las 
funciones para las cuales nosotros sabemos calcular las derivadas, 
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resulta fácil hallar la derivada de f(x). Con este fin se emplean 
las siguientes reglas que se demuestran (así como también las 
fórmulas 1—10) en el curso de las matemámicas superiores. 

1. La derivada de la suma de dos funciones es igual a la suma 
de las derivadas de estas funciones, es decir, 


AA + AY =41 00 +4209- 


2. El factor constante puede ser sacado fuera del signo de la 
derivada 
Laf = af). 


3. La derivada del producto de dos funciones se calcula por 
medio de la fórmula 


Da tos. CF =S O +A) 
4. La derivada de una fracción se calcula por la fórmula 


ES Į a LA -AOS 
1.6) PAE i 


La regla para calcular derivadas de un polinomio, expuesta en 
el $ 13, es un corolario de las reglas 1 y 2, y de la fórmula 2 
de la lista de las derivadas. 

Ejemplo 1. Hállese la derivada de la fracción 


3x = x+1 
To= F 


Haciendo uso de la regla 4, encontramos 


(Bx? — x + 17 (2x? + 5) — Qi? — x + D 2x? + SY 


Pr ES 


A continuación, según la regla para derivar polinomios 


(Bx? —x+1) =6x-1 


(2x3 + 5) =6x, 
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y. por esta razon, 


aa _ (xa? + S)(6x — 1) — (3x? — x + 1)6x? 
PIS AFA È 


_—6x% + 4x? + 6x? + 30x — $ 
aTM it 
(2x7 + 5) 


Ejemplo 2. Hállese la derivada de la función 
1 1 
Js= parsen g= +) 


Resolución, Según las fórmulas 2 y 9 y las reglas 1 y 2, 
tenemos i 


v 1 3 if-2 3 1 
AT ea TIMES 
Ejemplo 3. Hállese la derivada de la función 
fix) = 10*sen2x. 
Según la regla 3 y las fórmulas 3 y 4, tenemos 
fix) =(10'Y sen2x + 10*(sen2x) = 
= 10"sen2x In 10 + 10%-20052x = 
= 10* (sen 2x tn 10 + 2cos 2x). 


Las reglas indicadas permiten hallar las derivadas en muchos 
casos. Existe una regla más que es muy importante: la regla para 
calcular la derivada de una función compuesta. Esta regla se enuncia 
asi: 

Si la función y= fix) puede ser escrita en la forma y=Fl2). 
donde == p(x). entonces su derirada se define por medio de la 
fórmula 
dá Sx) = Fg). (17) 
donde == Q(x). 

Ejemplo. Hállese la derivada de la función y =sen(x3). Esta 
función puede ser escrita en la forma y =sen z, donde 22x. 
La derivada de la función F(:)=sen=: es igual a F'(z) = cosz, 
y la derivada de la unción ọ(x)=x° es igual a {x)= 3x? 
Aplicando la fórmula (77). obtenemos 

[sen x] = F()9'(x) = cos2-3x?. 


a 


Sustituyendo, en lugar de z, el valor z = x?, tenemos 
[sen(x)]/ = 3x7 cos(x?) 
La exposición más detallada del concepto de la derivada el 
lector puede encontrar, por ejemplo, en el libro “Matemáticas 


superiores para los primerizos” por Ya. B. Zeldovich, “Fizmat- 
guiz”, 1960. 


$19 
ELECCIÓN DE LAS PRIMERAS 
APROXIMACIONES 


Examinemos, ahora, el problema de la elección de las aproxi- 
maciones iniciales. Al resolver la ecuación f(x) =0, la búsqueda de 
la primera aproximación puede efectuarse por el método gráfico. 


a 
a) b) 19) 
Fig. 16 


Para ello es menester dibujar el gráfico de la función y = f(x) y encon- 
trar los puntos de intersección de este gráfico con el eje Ox 
(en estos puntos y=0, y, por tanto, f(x) = 0). 

Si, debido a ciertas circunstancias, resulta incómodo dibujar 
el gráfico (digamos, si la ecuación se resuelve con ayuda de una 
máquina calculadora), recurren al otro método de determinar la 
primera aproximación. Con este fin, se calcula el valor de la 
función para ciertos valores del argumento (por ejemplo, para los 
valores enteros del argumento ubicados dentro de los límites deter- 
minados). Si la función y = f(x) es continua (es decir, su gráfico 
no tiene puntos de discontinuidad), entonces, entre los valores a y b, 
en los que la función tiene signos contrarios (fig. 16, a), se encuentra 
la raíz de la ecuación f(x) =0 (si el gráfico de la función sufre 
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discontinuidades, puede ocurrir que ella pasa de salto de los valores 
negativos a los positivos, sin reducirse a cero (fig. 16, b)). Los valores 
a y b se pueden tomar por las primeras aproximaciones para la 
raíz de la ecuación f(x) = 0. 


Fig. 17 


Observemos que en este caso se puede omitir algunas raices 
de la ecuación. Así por ejemplo, la fig. 16, c ilustra el caso en que 
la función y = f(x) tiene en dos puntos los valores de un mismo 
signo, pero se anula entre estos valores. 

Hemos obtenido, de esta manera, dos puntos: a y b. Con 
el objeto de aclarar, cuál de ellos puede tomarse por la aproxima- 
ción inicial x en el método de Newton, examinemos la fig. 17. 
Las figuras 17, a y 17, b nos muestran que si la curva gira su 
concavidad hacia la parte positiva del eje de ordenadas, por la apro- 
ximación inicial se tomará aquel de los puntos a y b, en el que la 
función f(x) es positiva. La otra elección de la primera aproxima- 
ción puede incluso conducir a que el punto xz se encuentre fuera 
del segmento [a, b]. Del mismo modo. si la curva gira su concavidad 
hacia la parte negativa del eje de ordenadas, por la aproxima- 
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ción inicial se tomará el punto, donde la función f(x) es negativa 
(ñg. 17, c, d). 

Esta regla es cómoda para usar, si se conoce el gráfico de la 
función y = f(x). Si, en cambio, el gráfico de la función no está 
dibujado, se necesitan cálculos adicionales para determinar la 
dirección en que la curva gira su concavidad. Para ello se debe 
encontrar la segunda derivada de la función f(x). Se denomina 
segunda derivada de la función f(x) la derivada de la primera 
derivada. 

Por ejemplo, sea dada la función 


f0=x -4+3x-1, 
entonces, su primera derivada es 
f()=3x-8x+3, 
y 3u segunda derivada es 
f(x) =6x-8. 


En el curso de las matemáticas superires se demuestra que si en 
el segmento [a, b] la segunda derivada es positiva, la curva en 
dicho segmento gira su concavidad hacia la parte positiva del eje 
de ordenadas, Sí, en cambio, la segunda derivada en el segmento 
[a, b] es negativa, la curva gira su concavidad hacia la parte nega- 
tiva del eje de ordenadas. Al hacer uso de esta circunstancia, 
obtenemos la siguiente regla para aplicar el método de Newton: 

Supongamos que erm los puntos a y h la función f(x) tiene 
signos contrarios, con la particularidad de que en el segmento [a, b] 
la segunda derivada de la función f(x) es positiva. Entonces, por la 
primera aproximación x, se debe elegir aquel de los puntos a y b, 
en el que la función f(x) toma un valor positivo. Si en el segmento 
[a, b) la segunda derivada es negativa, por la aproximación inicial se 
escogerá el punto, donde la función f(x) toma un valo? negativo. 


; $20 
METODO COMBINADO PARA 
RESOLVER LAS ECUACIONES 


A) resolver las ecuaciones, con frecuencia recurren a la combi- 
nación de los métodos de cuerdas y de Newton. Si el gráfico de la 
función y=f(x) gira su concavidad hacia la parte positiva del 
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eje de ordenadas, los puntos a, y x, se hallan por medio de las 
fórmulas 


b~a 
ON ds 
fib) 
== 7 (m) 


Si el gráfico de la función gira su concavidad hacia la parte 
negativa del eje de ordenadas, el punto a, se halla por la fórmula 
(78), mientras que el punto x, se obtiene de la expresión 


me LA (80) 


Fig. 18 


Como lo muestran las figuras 18, a y 18, b, la raiz E de la 
ecuación f(x) = 0 se ubica, comúnmente, entre los puntos obtenidos 
a, y x. Aplicando, nuevamente, a estos puntos las fórmulas de los 
métodos de cuerdas y de Newton, se obtiene otro par de puntos, 
az Y Xn Cto. 

De este eS se obtienen dos sucesiones de puntos aj, az 

a Moo Y Xis Xa ce Xm -». que van aproximándose hacia 
la raíz buscada E desde las direcciones opuestas. La ventaja del 
método descrito consiste en lo que en este caso se obtienen 
los valores aproximados tanto por exceso como por defecto. 


Ejemplo. Resuélvase por el método combinado la ecuación 
x—senx-—05=0 


con un error inferior a 0,001. 
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Compongamos una tabla de valores de una funcion continua 
fix) =x —senx — 0,5. 


[x er o i 2 
Lu 0.659 -05 —0341 0,591 


De esta tabla se desprende que la raiz de la ecuación se ubica 
entre | y 2. Haciendo uso de las fórmulas 2 y 4 del $ 18, tenemos 


SA) = 1 cosx. 


Por ello, la fórmula de Newton toma cn nuestro caso la forma 


(81) 


Para determinar. cuál de los valores, 1 ó 2, hay que tomar 
por Xp, hallemos la segunda derivada de la función f(x). Según la 
fórmula 5. § 18. la derivada tiene la forma f (x/=senx. Mas, en el 
segmento [1, 2] la función sen x es positiva*. Por consiguiente, 
rigiéndonos por la regla citada anteriormente, por xo sc debe tomar 
el valor 2, para el cual también es positiva la función f(x) 

Según la fórmula (81) tenemos 

E 2-sen2-05__ 2-0909—0,5 
q 10082 7 140416 
Por otra parte, según la fórmula (78) tenemos 
2-1 
== —, m =1.366 
0,591 — (0.341) 

Siguiendo a aplicar las fórmulas ($1) y (78) al segmento [ar 

x]. obtenemos 


= 1,583 


a, = l — (— 0,341) 


1,583 — 1,000 — 0,5 


22 1.588 T 


1,583 — 
az= 13664011359 


A continuación, encontramos 
x, = 1,498. 
a, = 1,498. 
* Sen x es positivo en el segmento |V: x] = |0; 3141.. J 
Por lo tanto. sen v es también positivo en la parte [1. 2] de este segmento. 
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Asi pues, la raiz de nuestra ecuación con un error inferior a 
0,001 es 1,498. 


$21 
CRITERIO DE LA CONVERGENCIA 
DEL PROCESO DE JTERACIONES 


Emplearemos, ahora, el concepto de la derivada para obtener 
un nuevo criterio de la convergencia del proceso de iteraciones. 
Nos bará falta, para ello, una fórmula, llamada fórmula de Lagrange 
(Lagrange es un matemático del siglo XVIII, que vivió en Francia). 


Fig. 19 


Examinemos en el segmento [a, b] la curva y = f(x). Designemos 
con M el punto inicial de esta curva, y con N, su punto final, 
y tracemos la cuerda MN. El coeficiente angular de esta cuerda 
tiene por expresión 

kama tay EN 
cuerda = Ey = Gp 


ífig 19) Pero, MP=b—a, y PN = fib) — fía). por lo que 
k = 2h -sia 
cuerda a 
Destgnemos por Tun punto del arco MN que sea más alejado 
de la cuerda MN. Si por este punto trazamos una recta, parale- 


El 


lamente a la cuerda, ella será tangente a la curva; en el caso 
de intersección, en Ja curva habrían puntos más alejados de la cuerda 
MN, que el punto T. En otras palabras, la tangente a la curva, 
trazada en el punto T, es paralela a la cuerda MN y por esta 
razón tiene el mismo coeficiente angular que la cuerda. Más, el 
coeficiente angular de la tangente es igual a f(c}, donde c es la 
abscisa del punto T. Por ello, resulta válida la fórmula 


vi fib -sa 
A E 


Esta es la fórmula de Lagrange. Indiquemos que el punto c 
en la fórmula de Lagrange siempre se ubica entre los puntos a y b. 
La fórmula de Lagrange puede ser escrita también en la forma 


siguiente 
fib) -fid = f'(b — a). (83) 


(82) 


Volvamos ahora a la resolución de la ecuación x = (x) por 
el método de iteraciones. Supongamos que la aplicación y = p(x) 
transforma el segmento [a, b] en sí mismo, con la particularidad 
de que en este segmento se cumple la desigualdad leo] <q 
donde q es un número menor que la unidad, q < 1. En el segmento 
[a, b] tomemos dos cualesquiera puntos x, y x}. En este caso, 
los puntos p(x,) y Ọ (x2) también pertenecen al segmento a, b] 
Según la fórmula de Lagrange tenemos 


9(x2) — p(x1) = p(x — xi) 


donde c es un punto que se ubica entre x, y Xz, y, por tanto, 
ertenece al segmento [a, b]. Es válida la desigualdad 
9'(c)| <q <1, y, por esto, 

loxa- p) < q1x2 xil (84) 


La desigualdad (84) muestra que p(x) prefija una aplicación 
contraída. Pero, ya sabemos que si x— ẹ(x) es la aplicación que 
contrae el segmento [a, b] en sí mismo, para todo punto xy de este 
segmento la sucesión xo, Xp .»<» Xm ===» donde x+: = (Xp), 
converge la raíz de la ecuación x= p(x). De este modo, hemos 
demostrado el siguiente teorema: 

Teorema. Supongamos que la función y =p(x) define la apli- 
cación del segmento [a, b] en sí mismo y que en este segmento se 
cumple la desigualdad |p'(x)] < q, donde q < 1. Entonces, sea cual 
fuese el punto xp del segmento [a, b], la sucesión de puntos Xo, 
Xis -+e Xm <<<» donde Xp», =9(x,) converge hacia la raiz de la 
ecuación x = Q(x). 
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En términos generales el sentido del teorema demostrado consiste 
en lo que el proceso de aproximaciones sucesivas permite hallar 
aquellas raices E de la ecuación x = ọ(x), en las cuales se cumple 
la desigualdad le®l < 1. Puede decirse que estos puntos atraen 
hacia sí una quebrada que representa geométricamente el proceso 
de iteraciones (véase el $ 8), y los puntos, donde lo'(E)| > G 
repelen de sí la quebrada citada. 

Si la desigualdad |¢' (x)| <q <1 se cumple en todo el éje 
numérico, el proceso de iteraciones converge con cualquier elección 
de la aproximación inicial xy (véase el $ 10). 

Ejemplo 1. ¿Puede emplearse el método de iteraciones para 
resolver la ecuación 

cosx + Senx_,, 
E 


En el caso dado 
cosx + senx 
0. LEE 


Por eso 
~senx + cosx 


ten 
Pero, |senx| < 1, |cosx] < 1, por lo que 
Įsenx| + |cosx} at 
4 r 
y, por tanto, el proceso de iteraciones es aplicable. 


, —sen x + cosx 
bl — 


Ejemplo 2. ¿Es aplicable el proceso de iteraciones para resolver 
la ecuación 
x=4-2%? (85) 


La raíz buscada se ubica en el segmento [1, 2], puesto que la 
función continua y = x — 4 + 2* en este segmento cambia de signo: 


1-44+2'<0 y 2-4+2*>0, 
Pero, en el caso dado tenemos 
9 (x)= —2*In2. 


Acotemos la expresión 2*In2 en el segmento [1, 2]. Si 


1<x<2, entonces 2 < 2* < 4, 
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y, por esto, 
2In2 < 2*In2 < 41n2, 
De las tablas de logaritmos naturales (cuya base es el número 


e = 2,78...) tenemos In2 = 0,69... Esto significa que en el segmento 
[1, 2] se cumple la desigualdad 


1,38... < 2*ln2 < 2,76... 
y el proceso de iteraciones no es aplicable. 


Para aplicar el proceso de iteraciones, transformemos la ecua- 
ción (85). Escribámosla en la forma 


Y=é-x 
y tomemos logaritmos en ambos miembros según la base 2. 
Entonces obtendremos 
x = log: (4 — x). 


En este caso 
NE E EE A 
ro) =- gyni 


y en el segmento [1, 2] está cumplida la desigualdad 
A 1 J 
<= is! 


(El lector mismo puede deducir con facilidad esta desigualdad). 
Por esto, escrita la ecuación en esta forma, el proceso de ite- 
ración converge. 


A AAA 
$22 
RAPIDEZ DE LA CONVERGENCIA 
DEL PROCESO DE ITERACIONES” 


¡_xXIAAm—+ 


Hagamos uso de la derivada de la función p(x) para estimar 
la rapidez de la convergencia de iteraciones al resolver la ecua- 
ción x= ọ(x). Nuestra tarea es evaluar la velocidad con la cual 
disminuyen los errores % =é — x, de valores aproximados Xp, =-.. 
Xw «.. de la raiz E. 


* Este párrafo puede ser omitido en la primera lectura. 
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Observemos que son válidas las igualdades E =p(E) y xp+1 = 
= Q (x,). De éstas se deduce 


Anti = É = Xp = PE) — (xa). 


Pero, según la fórmula de Lagrange tenemos 
9) Phx) = P (0 (E — xn) = O (Ca) o 
donde c, es un punto ubicado entre los puntos x, y E. Por 


esta causa 
Onti = O (Ca) A (86) 


De la igualdad (86) se deduce la afirmación siguiente: 

Supongamos que E es la raíz de la ecuación x= Q(x), se ubica 
en el segmento [a, b]. Si en este segmento se cumple la desigual- 
dad |p(x)| <q < 1, y la aproximación inicial x, también está 
elegida en el segmento [a, b], entonces para todo n se verifica la 


correlación 
lanti l< glo l. (87) 
Efectivamente, de la igualdad (86) tenemos 
a: = |g'(c)lia l. 


Pero, el punto c, se encuentra en el segmento [a, b] (fig. 20) 
y, por lo tanto, 


lo (eni<e 


De aquí se deduce que 
laz! < qla]. 
Del mismo modo obtenemos que 
lasi =|9'(c2)]l021<gla1<g as! 
y, en general, 
länsi} < gla. 


Con esto queda demostrada nuestra afirmación. 
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Puesto que, para 0 <q < 1, la sucesión de números q, EA 
q’, ... tiende a cero, también a+, tiende a cero cuando n crece. 
En otras palabras, hechas las suposiciones indicadas anteriormente, 
los números Xy, ..., X» -.. Se aproximan al número č, con la 
particularidad de que la diferencia |§ — x,,,1 decrece con mayor 
rapidez que |a; įg". 

Del modo igual se puede demostrar que si en el segmento 
[a, b] se cumple la desigualdad 

le 91> 1, 


el proceso de iteraciones diverge. 

El proceso de aproximaciones sucesivas converge con una rapidez 
particular, cuando en el punto E la derivada de la función p(x) 
se reduce a cero. En este caso, el valor de ọ'(x) tiende a cero 
a medida que nos aproximamos a E. 

Dado que 
ln+11=19(c)11 0.1, 


la convergencia del proceso se hace más rapida con la apro- 
ximación al punto E. 

Ya nos tropezamos con esta circunstancia al extraer raíces 
cuadradas por el método de iteraciones. Recordemos que en esta 
operación sustituimos la ecuación x?=a por la ecuación 

2 2 
E =+ 2. Pero, la función arbitraria p(x) = z Z a 


la expresión 
(x? +a 2x — (x? +a) (2x _ 
4x? T 


x= es igual a 


g(x) = 
2x-2x - (x? +0)2 _ *-a 
E 4x 2x7 


(véase la regla 4 del § 18 y la fórmula (55) dei § 13), y, por lo 
tanto, 


Va? — a 
i = =——=0, 
9 Va WF 


De este modo, en el punto x=a la derivada de la función 
p(x) se reduce a cero y esto provoca la aceleración de la conver- 


6-71 
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gencia_del proceso a medida que nos aproximamos al punto 
x=ya, 

El fenómeno de la aceleración del proceso con la aproxima- 
ción a la raíz de la ecuación se observa también al emplear el 
método de Newton (el método indicado de la extracción de raíces 
es un caso particular del método de Newton). Efectivamente, ya 
indicamos que el método de Newton está relacionado con la susti- 
tución de la ecuación f(x) =0 por la ecuación 


Su 


y lá resolición posterior de la última por el método de apro- 
ximaciones sucesivas. En este caso tenemos 


SO) 
o{x) > 
pero 
o LON LOLI 
e0=:-[43] LAO 
-1 - EOE -SO - S04 
"0, ro 


Ya que en el punto E se cumple la igualdad f(&) = 0, enton- 
ces p'(Ẹ)=0. Y esto, según lo mostrado, asegura la aceleración 
de la convergencia del proceso de aproximaciones a medida que 
nos acercamos al punto E. 


$23 
RESOLUCIÓN DE LOS SISTEMAS, 

DE ECUACIONES LINEALES POR EL MÉTODO 
DE APROXIMACIONES SUCESIVAS 


Hasta este momento resolvíamos las ecuaciones con una incóg- 
nita. Pasemos ahora a la resolución de los sistemas de ecuaciones. 
Iniciemos con la consideración de los sistemas de primer grado. 
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Sean dadas m ecuaciones del primer grado con m incógnitas * 
NET E + GimXm = bys 


a1X3 + 22X2 +. + zoo = ba, (88) 


Ami Xy + Om2X2 + -o Om = Di 


Tales sistemas intervienen en diferentes aplicaciones. Por ejemplo, 
cuando los geodestas nivelan los resultados de las mediciones de las 
grandes partes de la superficie terrestre, ellos tienen que, a veces, 
resolver los sistemas de ecuaciones que se componen de muchas 
centenas de ecuaciones, Los mismos sistemas de ecuaciones aparecen 
ante los ingenieros, al calcular las construcciones de vigas, y ante 
otros especialistas. 

La resolución de tales sistemas por los métodos habituales 
(por ejemplo, usando el método de eliminación de las incógnitas) 
suele ser con frecuencia extremadamente dificultosa. Más cómodo 
resulta el método de aproximaciones sucesivas. Al principio mostre- 
mos con un ejemplo, cómo se hace esto. 

Sea dado el sistema de ecuaciones 

10x, — 2x3 + x3=9, 
1+5x,-x3=8, 


4x, + 2x2 + 8x3 = 
Es necesario hallar las incógnitas X,, xz, xj, con un error 
inferior a 001. 


Despojamos x; de la primera ecuación, x} de la segunda y 
x, de la tercera ecuación. El sistema tomará la forma 


X, =09 + 0,2x, — 0,Lx3, 
xz = 1,6 — 0,2%, + 0,2xy, (89) 
Xy =4—0,5x, — 0,2x2. 


Tomemos cualesquiera valores de xı, xz X3 a titulo de las 
aproximaciones iniciales, por ejemplo, hagamos x/%=0, x/%=0, 


* En estas ecuaciones las incógnitas se designan por 
las letras x}, Xz -<., Xw Y los coeficientes con las letras ap Aqui, el 
primer signo indica el número de la ecuación y el segundo, el número 
de la incógnita. Por ejemplo, as, significa que esto es el coeficiente de la 
séptima incógnita en la cuarta ecuación. 
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x4” =0. Sustituyamos estos valores en los segundos miembros de 
la igualdad (89) y tomemos los valores de xı, xz, x3 que se 
obtienen por los siguientes valores aproximados. Obtenemos 


xi” =09; 
x49 =1,6; 
xP =4. 


Los valores hallados los sustituyamos de nuevo en los segundos 
miembros de la ecuación (89). Obtenemos las aproximaciones 


xi? =09 + 0,2: 1,6 — 0,1 -4 = 0,82, 
xi? = 1,6 — 0,2-0,9 + 0,2-4 = 2,22, 
xP = 4 — 0,5-0,9 — 0,25: 1,6 = 3,15. 


En general, una vez hallados los valores de x/”, x{™, xf, 
las siguientes aproximaciones se calculan según las fórmulas 


xf0=09 + 02x4” — 0,1x7, 
xfP =16-0,2x/" + 0,2x5", (90) 
xfD=4-0,5x/” — 0,25xf”. 


Los resultados de los cálculos se exponen en la Tabla 2. 
Tabla 2 


Vemos que en los límites de la precisión prefijada se cumplen 
las igualdades 
x= x, 149 = 110, xÍ9 =x50, (91) 


Haciendo en las igualdades (90)n = 5 y tomando en conside- 
ración las igualdades (91), obtenemos que en los límites de la 
precisión prefijada 


xi? 209 +02x/2 — 0,1 x49, 
x4? = 1,6 = 0,2 x4" + 0,2 x49, 
xi? m4 — 0,5 x9 + 0,2514 
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A _Á EA -> > — 2 A 2 KK 


(en realidad, estas igualdades se cumplen exactamente, pero esto 
no es esencial). De aquí se deduce que los números xi’ = 1,00; 
x4” = 2,00; x{® = 3,00 son (en los límites de la precisión prefi- 
jada) las soluciones del sistema dado de ecuaciones. 

De un modo sumamente igual se procede en el caso general » 
Sea dado el sistema de ecuaciones (88). Despojemos x, de la 
primera ecuación, x; de la segunda ecuación, etc. Entonces, el sis- 
tema (88) tomará la forma 


x= x _ —E xo 
as ar an 
a 

a. Lao... — y 
a n % 
Bn e 

Xp = Om eps amz x= Am, m-1 AEA (92) 
Omm Onm Can a, 


Sean xf", ..., x4} ciertas aproximaciones iniciales para las in- 
cognitas Xi, -s Xm 

Sustituyendo estos valores aproximados en los segundos miem- 
bros de las igualdades (92), obtendremos las segundas aproxima- 
ciones x/?, ..., x4? para las incógnitas buscadas: 


b, a a 
AV me E gg E, 
A ar ai 
b a a; 
iD m 2 RL ra aai 
a2 a22 a22 
paap bu TNT] 


Del mismo modo, una vez halladas las n-ésimas aproxima- 
ciones x,™, x2 0", ..., x5 para las incógnitas, el siguiente paso se 


» El texto ulterior hasta el fin de este párrafo puede 
ser omitido en la primera lectura. 


efectúa según las fórmulas 


b, a, a 
PS Es A A 

A LN 831 

b. a a; 
x= n 1 o a o, 

422 422 a22 


Estudiemos, ahora, el problema en el que se trata no de la 
resolución del sistema de ecuaciones lineales por el método de 
aproximaciones sucesivas, sino al revés, resolviendo un sistema 
de ecuaciones se debe hallar un estado límite de cierto proceso de 
aproximaciones. 

Se tienen tres cubos. El primer cubo contiene 12 litros de agua, 
los cubos segundo y tercero son vacios. Del primer cubo se trasiega 
la mitad del agua en el segundo; luego, la mitad del agua se trasiega 
del segundo cubo en el tercero y de éste último la mitad del agua 
se trasiega en el primer cubo. Este ciclo se repite 20 veces. Hace 
Jalta determinar (con un error inferior a 0,0001 1), ¿cuánto agua, 
después de esto, habrá en cada uno de los cubos? 

Está claro que aquí se trata de las aproximaciones sucesivas 
hacía cierta distribución límite del agua. Esta distribución límite se 
caracteriza por lo que ella no cambia después de realizarse uno 
de los ciclos de trasegar descrito más arriba. Si al principio 
del ciclo en el primer cubo están contenidos x litros de agua, 
en el segundo y litros y en el tercer cubo, 12 — x — y litros (el volu- 
men general del agua no se altera durante el proceso de trasegar), 
entonces el ciclo descrito se da por medio de la siguiente tabla: 


l-r cubo 2-do cubo 3-r cubo 


Estado inicial x y 12=x-y 

Después del primer trasiego 3 3 +y 12-x-y 

Xx X4Y Y 

Después del segundo trasiego a +3 12-33 
ė 2 Y Eez 
Después del tercer trasiego ei 34 
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Para que este ciclo de los trasiegos no cambie la cantidad 
del agua en cada cubo, se deben cumplir las siguientes ecuaciones 


x_y 
643 7 


Al resolverlas, encontramos que x = 6, y = 3. Esto significa que 
la distribución límite se caracteriza por lo que en el primer cubo 
se encuentran 6 litros del agua, el segundo contiene 3 litros y en el 
tercer cubo habrán 3 litros del agua. 

Aclaremos, a qué velocidad la distribución dada del agua se 
aproxima a la distribución límite. Supongamos que el primer cubo 
contiene a litros y el segundo cubo, b litros. Realizado el primer 
ciclo, en el primer cubo habrán 


b 
a=6+5-3 04) 
litros, en el segundo cubo 
a,b 
b=3+3 (95) 


litros. 
Designemos a — 6cona, b— 3conf;a, — 6 mediante &,, y b, — 3, 
mediante f3,. Entonces, de las igualdades (94) y (95) proviene que 


a—-6 


B =b -35 


4 


YEN 3 5 
-a4 +i) -a-t 


tfa B 5 3 
A) 


38 


Por esto, si |a| <£ y |B]< e, entonces 


13 
lal < -g7 E= 025 


1 
B21 < 37E = 034e. 


Esto quiere decir, que dos ciclos de los trasiegos disminuyen 
los errores « y B por lo menos en tres veces. Por ello, después 
de 20 ciclos el error disminuirá por lo menos en 3*” = 70000 veces. 
Por lo tanto, realizados 20 ciclos de los trasiegos, con un error 
inferior a 0,0001 1, en el primer cubo habrán 6 litros de agua, 
en el segundo cubo, 3 litros, y en el tercer cubo, 3 litros. 
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RESOLUCIÓN DE LOS SISTEMAS 
DE ECUACIONES NO LINEALES POR 
EL MÉTODO DE APROXIMACIONES 

SUCESIVAS 


El método de aproximaciones sucesivas (de iteraciones) se emplea 
también para la resolución de algunos sistemas de ecuaciones no 
lineales. 

Consideremos, por ejemplo, el sistema de ecuaciones 


s=2+ ŽAL, 

, 69) 
e x+y 
y =14+—3 


A título de las primeras aproximaciones tomemos xo =Q, 
Yo = 0. Sustituyendo estas aproximaciones en lugar de x e y en los 
segundos miembros de las ecuaciones, obtendremos las aproxima- 
ciones siguientes: x,=2, y, = Í. Sustituyendo estas aproxima- 
ciones en los segundos miembros de las ecuaciones (96), obtenemos 
PRL 
a” 2,25, 


241? 
20 


x=2+ 


y=l+ 


= 1,5. 
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Continuando este proceso, obtendremos 
225 +1,15 _ 


9.24 2% 2,31, 
Yy=1+ e 1,18, 
aA PA 233, 
ya=i A 1,18, 

x=2+ A L 233, 
ys=l+ PA, 1,18. 


Vemos que con un error inferior a 0,01 se cumplen las igual- 
dades xy = xs = 2,33 y ya = ys = 1,18. Por esto, con la precisión 
indicada tenemos x = 2,33 e y= 1,18. 

En general, si se da el sistema de ecuaciones 


ia 


(97 
y =y y), on 
donde p(x, y) y y (x, y) son ciertas funciones, entonces eligimos las 
aproximaciones iniciales xy e yo, los sustituimos en los segundos 
miembros de las ecuaciones (97) y luego hacemos los cálculos 
según las fórmulas 


Xara = Om Y) (98) 
ner = Ym Ya) + 

Si para cierto número n, con el grado dado de precisión se cumplen 

las igualdades X,+1 Y Xm )n+1 Y Xp Entonces con el mismo grado 

de precisión tenemos x 2 Xp y Y Yw 

De esta misma manera se resuelven los sistemas de ecuacio- 
nes con tres y mayor número de las incógnitas. 

Aclaremos, ahora, en que condiciones puede ser garantizada 
la convergencia del proceso de aproximaciones sucesivas al resolver 
los sistemas. 

Vamos a considerar que las funciones p(x,y) y v(x,y) en el 
sistema de ecuaciones (97) están definidas en cierto dominio ce- 


1-7 
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trado y limitado Ð en el plano x, y. En otras palabras, consi- 
deremos que el dominio D se encuentra enteramente dentro de un 
cuadrado y que le pertenecen todos sus puntos limites. Como ejemplo 
de tales dominios pueden servir un círculo, un polígono (tomado 
junto con una quebrada límite), una elipse, etc. Además, vamos a 
considerar que las funciones p(x, y) y Y(x, y) son continuas en el 
dominio D. Las funciones 9 (x, y) y Y{x, y) determinan la aplicación 
del dominio D en algún otro dominio ubicado en el mismo plano, 
Para hallar el punto en el cual se transforma el punto Mo (Xo, Yo) 
del dominio D, es necesario sustituir sus coordenadas, en lugar 
de x e y, en p(x, y) y Y(x, y). Los resultados de la sustitución 
nos darán las coordenadas de la imagen del punto Mọ. Por 
ejemplo, si 
oa =y, 


yix, y) =2xy, 


entonces el punto Mo (t, 3) se transforma en el punto Ne (10, 6). 

En lo sucesivo designaremos con la letra ® la aplicación dada 
por las funciones p(x, y) y Y (x, y) y mediante ®(M), la imagen 
del punto M para esta aplicación. Por W(D) designaremos la imagen 
de todo el dominio D. 


Er 


D 


Fig 21 


Supongamos que la aplicación ® transforma el dommio D 
en una parte de éste D, =D(D). Entonces, podemos realizar 
reiteradamente la misma aplicación. En este caso, el dominio D, 
se transformará en su parte D} = d(D,), la cual también se ubica, 
por supuesto, en el dominio D. Continuando este proceso, obten- 
dremos un sistema de los dominios DÐ, D,, ..., D,, ..., encajados 
uno en otro (fig. 21). 
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Llamaremos la aplicación D contraidu. si existe un numero 
q tal que 0<g<l y que para dos puntso cualesquiera M, y M2 
del dominio D se cumple la desigualdad 


rM). D(M))<qr(M,. Mal. 


Aquí con r(M, N) está designada la distancia entre los puntos 
My NM. 

Tgual que en el caso de una sola variable tambicn se demuestra 
ta siguiente afirmación. 

Supongamos que la aplicación Y transforma el dominio D en una 
parte de éste y que es contraida. Entonces, en el dominio D existe 
un mico punto N tal que N =D(N). Este punto pertenece a todos 
los dominios D, Las coordenadas €, € del punto N satisfacen el 
sistema de ecuaciones 197). es decir 


¿=09( n 
n= m. 


Como en el caso de una sola variable. los valores Ų y n se 
calculan aproximadamente por el método de iteraciones, Si 
Mo (Xas vn) es un punto cualquiera del domino D y Mn+, = D(M,) 
(es decir, Xuri = O(m Yad Yna 1 = WAX: Fm). Entonces. la sucesión 
de los puntos My. My. ~ M, ~ converge hacia el punto inmóvil 
N de la aplicacion ® 


25 
DISTANCIA MODIFICADA 


wn 


La condición de que la aplicación D cs contraida nos presta 
el criterio suficiente para la convergencia del proceso de iteraciones. 
Más, este criterio no es necesario. La aplicación @ puede no ser 
contraida y el proceso de iteraciones. no obstante, es convergente 
Por ejemplo. la aplicación O, dada por las funciones p(x. y) = 


= 44 2v, dx 1)=3 +. no es contraida, Si tomamos A (8, Oh. 
B (8, 4). entonces 


r(A. B)=4 D(4)=(L, 4) D(B)=(9,4) 


HO(A) 0(B))=8 >r(4. B} 
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Sin embargo, cualquiera que sea el punto Mo eligido, la sucesión 
de los puntos Mo, My, .... My ... converge hacia el punto 


E id 
nfs}, at), 

En ciertos casos se logra establecer la convergencia del proceso 
de iteraciones al modificar la manera en que se determina la distancia 
entre los puntos de un plano. El hecho es que la distancia entre 
los puntos puede medirse a la manera diferente. Para un viajero 
es natural medir la distancia en las unidades del tiempo que se tarda 
para llegar al punto B, partiendo del punto A. En este caso, 
la distancia entre los puntos A y B en la fig. 22, a se determina 


B 


Fig. 22. a, b 


por la suma de longitudes de los segmentos AC, CD y DB (para 
llegar al punto B, partiendo del punto A, es necesario alcanzar 
el puente CD, pasar este puente, y del punto D llegar hasta el 
punto B). Si el movimiento en un plano sólo puede realizarse 
en dos direcciónes mutuamente perpendiculares, entonces por la 
distancia entre los puntos A y B en la fig. 22, b se debe tomar 
la suma de los segmentos AC y CB. En otros casos es más có- 
modo considerar como la “distancia” entre los puntos A y B la longitud 
del mayor de los segmentos AC y CB. Pueden también idearse 
otras definiciones de la “distancia” entre los puntos de un plano. 
(Los mayores detalles sobre diferentes modos de medir las distancias 
se puede sacar del libro “¿Qué es la distancia?”, escrito por 
Yu, A. Schreider.) 
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Habitualmente se requiere que la distancia r(4, B) entre los 
puntos posea las siguientes propiedades: 

1. La distancia r(A, B) entre dos cualesquiera puntos A y B 
es no negativa, con la particularidad de que ella es nula sólo 
en el caso cuando A y B coinciden, 

-2. Para dos cualesquiera puntos A y B se cumple la condición 


de tir 
le simetría r(A, B)=r(B, A). 


3. Para tres cualesquiera puntos A, B, C se cumple la desigual- 
dad triangular 
r(A, B} < r(A, C)+r(C, B). 


Si en algún conjunto de objetos queda determinada una distancia 
provista de estas propiedades, dicho conjunto se denomina espacio 
métrico, y los elementos del conjunto son puntos de este espacio. 
El papel de los puntos de un espacio métrico incluso pueden 
desempeñar las funciones. Para las funciones continuas p(x) y y (y) 
en el segmento [a, b] la distancia se determina como el valor 
máximo en este segmento de la función | y (x) — Y (x)|: 


rie. y) = máx 1900) — (01 


Ya vimos más arriba, que un plano puede ser transformado 
en un espacio métrico por los métodos diferentes: para todos 
los métodos de determinar las distancias, citados anteriormente, 
se cumplen las condiciones 1) — 3). 

He aquí un ejemplo interesante en que las condiciones 1) y 3) están 
cumplidas, mientsas que la condición de simetria 2) no se cumple. Vamos 
a medir la distancia entre los puntos A y B en un terreno montañoso 
en las unidades de tiempo necesario para llegar del punto A al punto B. 
Puesto que el tiempo para subir la montaña se diferencia del tiempo 
del descenso, entonces, r(A, B) + r(B, A). 


Resulta que para la convergencia del proceso de aproxima- 
ciones sucesivas al resolver el sistema de ecuaciones 


ná } (99) 
y= Vx y) 
es suficiente que se cumpla la siguiente condición: 
La aplicación ®, prefijada por las funciones p(x, y) y W(x, y), 
transforma el dominio D en si mismo y es contraída respecto 
a siquiera una sola “distancia” r(A, B). En otras palabras, debe 
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q€_ - > === o —_—2—<—Á 


existir tal número q que 0 < q < 1, y para cualesquiera dos puntos 
Mı y Ma de D se cumple la desigualdad 


r(d(M y), DM) <gr(M;, Ma). 
“Tomemos, por ejemplo, las funciones p(x, y)=1+2y y 


y(x y) =3 +5 Resulta que la aplicación D, definida por ellas, 


es contraída y tiene el coeficiente 1/2, si la distancia r(4, B) 
entre los puntos A(X1, yı) Y B(xz, Ya) es determinada por medio 
de la fórmula 


1 
r(A, B)=|50% = xı) + 20 -v A 


z — x1) — 2% — ya). 


Por ejemplo, para los puntos A(8, 0) y B(8, 4) tenemos 
r(A, B)= 16, y para sus imágenes D(A) y ®(B) tenemos 


r(D(A), D(B)) = 8. 
Realizada la aplicación D, la “distancia” se disminuyó dos 


veces. Precisamente por esto se explica que el proceso de itera- 
ciones para la resolución del sistema de ecuaciones 


x=1+2y, 
x 
y=3+3 


converge, aunque la aplicación ® no es contraida respecto a la 
distancia ordinaria (véase la página 91). 


$ 26 
CRITERIOS DE LA CONVERGENCIA 
DEL PROCESO DE APROXIMACIONES 
SUCESIVAS PARA LOS SISTEMAS 
DE ECUACIONES LINEALES 


Apliquemos el criterio de convergencia obtenido a los sistemas 
de ecuaciones lineales. Al elegir de manera diferente las “distan- 
cias”, obtendremos los criterios de convergencia para estos sistemas 
expresados mediante las propiedades de sus coeficientes. 


95 


Examinemos primero un sistema de dos ecuaciones lineales 
con dos incógnitas 
x+ ay =bn 
a,x + ay = bi: } uoo 


az1X + a22 = bz- 


Sea a,1 40 y 7, # 0. Resolvamos la primera ecuación respecto 
a x, y la segunda, respecto a y. Obtendremos el sistema de ecua- 
ciones 


b a 
a 
as ar 
b: a 
ME 
da 2 
b a b, 
Para abreviar, hagamos —- =P -75 s0, —2= Ba 
ar 4h a2 
a, 
-2s 
4 


Entonces, el sistema tomará la forma 


x= y+ Br (100) 
y =x + Pa 
En este caso las funciones que prefijan la aplicación ®, tienen por 
expresión 
oix, y) =a +B Yx y) = 22x + Pz. 
Actaremos, cuáles deben ser los coeficientes a, y % para que 
esta aplicación sea contraída. 


Como sabemos, la distancia r(A, B) entre los puntos 
Ax. ya) y Blxa, ya) se expresa por medio de la fórmula 


A A a 
HA, B)= Vix = x) +277. 


La aplicación ® pasa el pumo A en el, punto 4, (4y, + Bo 
aax, + Ba) y el punto B,.en el punto B,(2,y2 + Bi %2x2 + Ba). 
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La distancia entre estos puntos se calcula según la fórmula 


riAy Bi) = (y2 = ayı)? + (097 = 02x1} = 
= Gb. ~ y + ax —x Y". (101) 


El mayor de los números |a,| y |æ] lo designemos con q: 
q = máx (jail, 1221). 
Entonces, de la fórmula (101) se desprende la desigualdad 


r(Ar, Bj) < q Vba = x1}? + (2 ~ y1}? =gr(4, B). 


Quiere decir que si q < 1, la aplicación ® es contraida en todo 
el plano. Entonces, como sabemos, el proceso de aproximaciones 
sucesivas converge. 

Hemos, pues, demostrado que si 


måx(|æ l, (a21) < 1, (102) 


el proceso de aproximaciones sucesivas, al resolver el sistema de 
ecuaciones (100) siempre converge. 

Expresemos el criterio de convergencia obtenido directamente 
en términos de los coeficientes del sistema de” ecuaciones (100). 
Para esto recordemos que 


Sustituyendo estas expresiones en la condición (102), obtenemos la 
siguiente deducción: 

Para que el proceso de aproximaciones sucesivas, al resolver el 
sistema de ecuaciones lineales 


ayx+ajyy=b,, 
azıx + a22Y = bz, 


sea convergente, es suficiente que se cumpla la condición 


4 a a: 
máx 22, CA 
aii a22 


Esta condición implica que los coeficientes diagonales deben ser 
mayores que los coeficientes no diagonales ‘dispuestos en la misma 
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línea. Por ello, por ejemplo, al resolver el sistema 
x-3Jy=-11, 
6x+y=10 
la primera ecuación tiene que ser resuelta respecto a y, y la segunda, 
respecto a x: 
Ml. 2d 


i C 
a a eegen 


A veces resulta útil transformar de antemano el sistema de 
ecuaciones, habiendo sustituido las incógnitas x e y por unas 
incógnitas proporcionales. Tomemos, por ejemplo, el sistema de 
ecuaciones 

2 = 14. 
12x+ y= 14 } (103) 


3x-2y=-1. 


Para este sistema 


812 az 131.3 
max ası |’ | a22 ) + max 2” ij- 


y, por tanto, no se cumple el criterio suficiente par la conver- 


gencia del proceso de aproximaciones. Pero, si hacemos xaz 


entonces obtendremos el sistema 
4z+y= l4, } 


z-2y= -1, 
para el cual 


e(l) 3) 


Esto significa que el sistema (104) puede resolverse por el método 
de aproximaciones sucesivas. 

Por supuesto, la idea de emplear el método de aproximacio- 
nes sucesivas para resolver tales sistemas sencillos como es el 
sistema (104) no va a ocurrir a nadie. Sin embargo, este método 
suele ser muy útil para los sistemas que tienen gran número de 
incógnitas. Los criterios suficientes de la convergencia para el 


(104) 


A12 
âi 


la 
la 
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sistema 
E +... + dinXy = bi, 


21X1 + 032X232 +... danX, = bz, (105) 


Da 


GX HO X2 +... + aX 


se establecen de una manera aproximadamente igual que para 
el sistema de ecuaciones con dos incógnitas. No obstante, en este 
caso los diferentes métodos de hallar la distancia entre los puntos 
conducen a los resultados diferentes. A saber, subsiste la siguiente 
afirmación. 

El proceso de aproximaciones sucesivas para el sistema (105) 
converge, si se cumple una de las siguientes condiciones: 


(106) 


(107) 


(108) 


njat 


Una raya en las sumas (106) y (107) significa que está omi- 
tido el sumando para el cual i æj. En la suma (108) el signo (k) 
significa que están omitidos no sólo los sumandos para los cua- 
les i = j, sino también aquellos para los cuales i = k. Observemos 
que la condicion 3) se cumple a ciencia cierta, si 


n 


y 


JEL 


2 


SE <a, (109) 


Qu 


donde la raya signífica que están omitidos los sumandos para los 
cuales ¿=¡j En la mayoria de los libros sobre las matemáticas 
de cálculo esta condición interviene precisamente en la forma (109). 
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Como ya hemos indicado, las condiciones (106) — (108) se obtie- 
nen si exigimos que la aplicación 


sea contraída respecto a cierta distancia. 
A saber, la condición (106) corresponde a la distancia 


r(A, B) =máx(|x, = Yal, -< [Xn — Yal), 


la condición (107) corresponde a la distancia r(4,B)= 
n 


= Dlx — yıl, y la condición (108) corresponde a la distancia 
ra 


r(4,B)= |/ „)? entre los puntos A (X1, -- , Xy) Y BOn +... 


Yn)- 


Como en el caso de dos variables, a veces resulta util sustituir las 
incógnitas x + X, por las incógmtas proporcionales a éstas, y, = p1X,, 
+... Yn = Pnn donde p,>0,. ., Pa > 0. 

En este caso, las condiciones (106), (107), (108) tomarán la forma 


n mY Per, (106) 
; P 
na 
2) max X aA ea, 107) 
i an | Py 


q k 
3) max y tu (1081) 


i jea 
En particular, cuando p; = |a4], estas condiciones toman la 
forma 


1) máx EA <t, (106) 
i ajy 


ja 


2 2u m 
) mix Y a e (107) 
(=1 
5) ma Y "| 
máx E <1 (108°) 


i J= 


A titulo de ejemplo examinemos el sistema de ecuaciones 
x — 06y — 0,5z = — 26, 
—0,2x + y -0,42 = 3, 
—0,1x + 0,5y + 2 = 3,9. 


Las condiciones 1), 2), como también 1), 2’) para este sistema 
no se cumplen. Del modo igual para él no tiene lugar la desi- 
gualdad (109): la suma de los cuadrados de elementos no diago- 
nales es igual a 1,07. Pero 


K) 
máx 
k ; 


y, por tanto, el sistema puede ser resuelto por el método de 
aproximaciones sucesivas. 

Indiquemos que todas las condiciónes, enunciadas más arriba, 
son suficientes para la convergencia del proceso de aproximaciones 
sucesivas, pero de ninguna manera no son necesarias. Eligiendo 
otros métodos para medir la distancia entre los puntos y apuntando 
la condición de contracción, obtenemos otras condiciones de la 
convergencia. Sin embargo, ahora no vamos a considerar este 
problema. 

Para los sistemas de ecuaciones lineales son válidas las mismas 
observaciones que hemos mencionado en el § 5. Por ejemplo, 
el resultado de las aproximaciones no depende de la elección de 
la aproximación inicial. Por esto, si en el transcurso de los cálculos 
fue cometido un error, éste no desvaloriza los cálculos ulteriores 
sino sólo demora la obtención del resultado final. 

Existen varios métodos de las aproximaciones sucesivas para 


ay)? 


uj = máx (0,6? + 0,5? + 0,2? + 0,42; 0,6? + 0,5? + 
gi 


+0,1? + 0,5%; 0,22 + 0,4? + 0,1? + 0,52) = 0,87 < 1, 
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resolver los sistemas de ecuaciones lineales. Por ejemplo, en algunos 
de ellos, al hallar el valor aproximado x¿"**, para encontrar xf"*", 
sustituyen los valores x"*9, xi, xf”, ..., x4; a continuación, 
para encontrar xj”*' sustituyen los valores x/*0, xP+D, x49, 
an Xi), etc, La descripción de todos los métodos posibles de 
las aproximaciones para los sistemas de ecuaciones lineales podría 
constituir el objeto de otro libro. 


$ 27 
APROXIMACIONES SUCESIVAS 
EN LA GEOMETRÍA 


Hemos descrito las aplicaciones del método de aproximaciones 
sucesivas para resolver las ecuaciones y los sistemas de ecuaciones. 
Este método también se emplea en algunas ramas de la geometría, 
por ejemplo, al calcular la longitud de una circunferencia. Como 
sabemos, para calcular la longitud de una circunferencia, primero 
se halla el perímetro del cuadrado, inscrito en ésta, luego el 
perímetro del octógono regular inscrito, luego el perímetro de un 
polígono regular de 16 lados, etc. El límite de estos perímetros es 
precisamente la longitud de la circunferencia, Con ello, cada peri- 
metro ulterior se calcula con ayuda del perímetro anterior, Esto se 
hace del modo siguiente: 

Designemos por A, un lado del polígono regular de 2" lados 
y por P,, el perímetro de éste. Por ejemplo, A, es un lado de un 
cuadrado, y por lo tanto, A2 =RY2, Pi = 4R V2. Supongamos 
que ya encontramos P,. Entonces, evidentemente, 

P 


An= > 


En la geometría se demuestra que el lado az, de un polígono 
regular inscrito de 2n lados se expresa a través del lado a, del 
poligono regular inscrito de n lados y del radio R de la circun- 
ferencia según la fórmula » 


* Esta fórmula se deduce con mayor facilidad con 
ayuda de la trigonometria, Es evidente que si a, es el lado de un polígono 
regular inscrito de n lados, y a es el fado de un poligono regular 


al 
am=R |/2- |/a- 5. (110) 


Esto quiere decir que el lado A+, del poligono regular de 
2"+! lados se expresa mediante el lado A, del polígono regular 
de 2” lados según la fórmula 


Ansi =R -1-3 


B; 
Como A, = Fry A 


Peth , entonces 
2 
Prs=21R |/2- |/4- m (110) 
La sucesión de los números Pz, Py, ..., Pa ... tiende a la 


longitud de ła circunferencia, es decir, al valor 2HR. Por esta 


inscrito de 2n lados, entonces 


r R 
a = 2Rsen Y 4, = 2R sen 3, 


(véase el dibujo). Puesto que 


x cosa 
senj = Ñ 


ES 
a pe 
2R += E 
sen 5 
rd a} 
=R a ART R |/2- |/a A 
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razón, la fórmula (110) puede considerarse como la fórmula para 
el cálculo del valor de 2R por el método de aproximaciones 
sucesivas. Empleando este método, se puede hallar el valor de n 
con cualquier número de signos decimales. 

Existe otro método del cálculo aproximado de n, denominado 
método de perímetros iguales. Al aplicar este método, el poligono 
regular de 2" lados se sustituye por el polígono regular de 2"*' lados 
que tiene el mismo perímetro. Designemos por Í la apotema del 
polígono de 2” lados y por rm el radio de la circunferencia circuns- 
crita. La apotema del poligono regular de 2"*' lados que tiene el 
mismo perímetro que el poligono de 2" lados la designemos por 
1,+1 Y el radio de la circunferencia circunscrita, por r,+, 

Sea AB (fig. 23) un lado del polígono de 2” lados, inscrito en la 
circunferencia de radio r, Unamos el centro C del arco AB con 


£ 


CAS, 


0 
Fig. 23 


los puntos A y B y tracemos la linea media DE del triángulo ACB 
Es evidente que el ángulo DOE es igual a la mitad del angulo 40B 
qm. 


Por esto, DE es el lado del poligono regular de 2"*' lados inscrito en 


[i 
5 4B. entonces el 


la circunferencia de radio OD. Puesto que DE = 


perímetro del poligono de 2"** lados es igual al perimetro del 
poligono de 2° lados. Esto quiere decir que r,+, = OD. l,.¡ = OK 
Es fácil calcular que 


l: =0K = (11) 
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Luego, del triángulo rectangular ODC encontramos 
fażi = folati- (112) 


Las fórmulas (111) y (112) expresan r„+1 y lh+1 en términos de 
EA 

Al crecer n, los perímetros de los polígonos no cambian y los 
números r, y [, van aproximándose a un mismo límite. Este límite 
es igual al radio de la circunferencia cuya longitud es igual al perí- 
metro de nuestros polígonos. Si elegimos el polígono inicial de una 
manera tal que su perimetro sea igual a dos, entonces tanto r, 


como |, tienden al número A 
1 


E nA 
imei lahs 


Por ejemplo, si tomamos por el primer poligono un cuadrado 


2 
de lado 1/2, entonces rı = i lh =} Por esta causa subsiste la 


2 
siguiente afirmación: si hacemos r; = VEN h= A y calculamos los 
valores raer h41, n=2, 3, -.. según las fórmulas (111) y (112), 
entonces 

i í 1 

lim r, = lim h = —. 

1 nx El 

Por medio de estas fórmulas podemos hallar aproximadamente 

el valor de Ł Con este fin se debe realizar los cálculos hasta 
que los valores de r, y l, coincidan en los límites de la preci- 
sión escogida, Este valor general de r, y l, será el valor de ; con 
la precisión prefijada. 


$2, 
CONCLUSION 


Hemos conocido en este libro el empleo del método de aproxi- 
maciones sucesivas en la resolución de diferentes problemas: la ela- 
boración de los planes, la extracción de raices, resolución de las 
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ecuaciones, cálculo de la longitud de la circunferencia. Con estos 
problemas no se agota la diversidad del uso de este método. 
Muchos problemas llevan a las ecuaciones diferenciales (en las que 
figuran las derivadas de las funciones buscadas), las ecuaciones 
integrales, y a las ecuaciones del tipo aún más complejo. El método 
de aproximaciones sucesivas también constituye uno de los métodos 
más acertados para la resolución aproximada de dichas ecuaciones, 
Por supuesto, en estos casos el uso del método de aproximaciones 
sucesivas es mucho más complejo en comparar con el caso de las 
ecuaciones algebraicas. No obstante, se puede afirmar que sin el 
empleo del método de aproximaciones sucesivas no sería posible 
resolver ni uno de aquellos grandiosos problemas fisicos y técnicos 
que hoy día se resuelven. Este método es de amplio uso en los 
cálculos, del movimiento de los sputniks, de los reactores nucleares, 
en la exploración de la estructura del átomo. Sin embargo, el 
relato sobre el empleo del método de aproximaciones sucesivas 
fuera de los límites de las matemáticas elementales no es la 
tarea de esta obra. 
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EJERCICIOS 


Para que el lector pueda comprobar en qué grado él se apoderó de 
los métodos de la resolución aproximada de las ecuaciones, expuestos en 
este libro, damos a continuación unos cuantos ejemplos de la resolución 


aproximada de ecuaciones. 


Resuélvanse por el método de iteraciones 


las ecuaciones * 
1 
L= s 
ay 8. x? = sen x, ; ' 
4. x n(x+ 1}. 
2. x= (x + 1) 9. x= arsen ŽEL, 

3 4 15. nx=4- x. 
taws PEL, 10. x = cosx. 16. Inx =2 ~ x. 
xi 1 A 
4. e 1l. x= 17. x2 =h 
DERE ya aosie 18. Inx =0,1 x 
S x=y5= x 12. x= 1+ <p senx 19. tgx = lgx. 


6. 4- x=tgx. 


LES 
T. x? = genx, 13. x= t gwta MAZA 


Resuélvanse por el método de Newton 
las ecuaciones. 


21. x-$x+1=0 24. 5 4+5+1=0 
22. x? ~ 9x? + 20x — 11 =0. 25. senx + x= l. 
23, xX -3è +11 =0. 26. x? — l0lgx — 3 = 0. 


27. Haciendo uso del método de las aproximaciones sucesivas, resuél» 
vanse con un error inferior a 0,001 los sistemas de ecuaciones 


0,2y — 0,1 z + 0,898, 
03x + 0,152 + 1,383, 
z =0,25x —04y + 3,677; 


x= 


ai. VX F 2y — 0,710, 


x = Esen(x + y) + 0.336, 
» f 
y = — sen(x — y) + 0.362; y=\y-x+!. 


* En unos de los ejemplos dados es necesario, pri- 
meramente, reducir la ecuación al tipo x= p(x). 


107 


RESOLUCIONES 
1. Hagamos px) = pr Entonces vi) = nr Tenemos 


1 
$(0)=1>0, p(1) O 1. Por ello, el segmento [0, 1] contiene la raíz 


de la'ecuación dada. Sin embargo, no podemos aplicar el método de 

aproximaciones sucesivas al segmento dado, puesto que |q'(0)|=2> 4 

Para hacer el segmento más estrecho, indiquemos que Q(0, 4)= 
4 


= qog > 04. Por esta razón. la ecuación tiene raiz en el segmento [04. 1). 


Además, |9'(x)| «Ír< 1, siempre que 04<x<1, y, por lo tanto, 


puede emplearse el método de aproximaciones sucesivas. Haciendo x; = 0,4, 
obtenemos, después de 11 pasos de la aproximación, que xi, = Q(x11) = 
= 0,4655, 

Por consiguiente, con un error inferior a 0,0001 tenemos x = 0,4655. 


2, Hagamos (x)=(x + 1). Entonces, g'(x) =3(x + 1)?, Tenemos 
0(-2) =-=1> -2,p(-3) = —8 < —3, Por lo tanto, el segmento [—3, —2] 
contiene la raiz de la ecuación dada. Sin embargo, no podemos aplicar 
el método de aproximaciones sucesivas, puesto que en el segmento [ 3, —2] 
tenemos |q'(x)| > 1. Escribamos la ecuación dada en la forma 

3 


x=|px- 1 


> 1 
Entonces, yip =|Y/X— 1 y Yx) = ——. En el segmento [-3, +2] 
3p 
t 
tenemos |y (x)| < ——< i. y, por lo tanto, puede emplearse cl método 
3px 


de aproximaciones sucesivas. Al hacer x, = —2, obtenemos xe = Y íx) = 
= —2,325, Esto quiere decir que con un error inferior a 0,00) tenemos 
x= - 2325. 


3. Hagamos p(x)=4+ |/ ==; - Tenemos 


ix + I 


La fig. 24 muestra que la recta v= x corta la curva y =4+ | Ei 
en dos puntos, ubicados, respectivamente, en los segmentos [ — L, 0] y [4,5] 
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En el segmento [4,5] tenemos | 9'(x)] < —< 1. Haciendo x,=4, 
15/45 


tenemos x, = Q(xs) = 4,870. Esto quiere decir que con un error inferior 
a 0,001 obtenemos x, = 4,870. 

En el segmento [—1, 0] el método de aproximaciones sucesivas no 
puede emplearse directamente. Escribamos en este segmento la ecuación 


A E STIR 
cn la forma (14) = 7» de donde (Xx +1 y x 
f Ea 
SRT! 
i x= Í dos A x 
Aqui, Y) = =$ -ly vo=- . Está claro que para —1 < 


<x<0 tenemos |y'(x)| <35 <l, y ahora podemos aplicar el método 


de aproximaciones sucesivas. Al hacer x, = 0, tenemos xz = Y (x2) = 0,9840, 
Por consiguiente, con un error inferior a 0,0001 tenemos x= 0,9840, 
Hemos encontrado dos raices: x = —0,9840, x = 4,870. 


Fig. 24 


.. 4. 

4. Aqui, 9, (1)=2+ | x y q2(x)=2— | x. De la figura 25 se ve que 
Es 

la ecuación x = 2 + |/x tiene raíz en el segmento [3, 4]. En este segmento 


19091 = <1 
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Empleamos el método de aproximaciones sucesivas Al hacer x; =4, 
tenemos x4 = p(x,) = 3,353. Por lo tanto, con un error inferior a 0,001 
tenemos x = 3,353, 


y 


yz 


Fig 25 


._ 
Ahora resolvamos la ecuación x=2-—]'x Tiene la raíz x= l. Así 
Pues, las raices de la ecuación son iguales a 1 y a 3,353. 
=$ 
5. Aqui, ọ {x)= %5 =x y 0()=—__——- 
35- x? 


Tenemos 3 š 
ol) =yV4<2 00)=13>1 
Por esto, en el segmento [1, 2] existe la raiz de la ecuación. En este 


segmento |9'(x)| < 


q Al hacer x, = t, tenemos xs = p(xs) = 1,516. 

31 

Por consiguiente, con un error inferior a 0,001 tenemos x = 1,516. 
6. Escribamos la ecuación en la forma 


x =arctg (4 — x). 


Aqui, p(x) = arctg (4 — x). Tenemos p(1) = arctg 3 = 1,25, 4 (2) = arctg2 x 

æ 1,10. De aqui se deduce que existe la raiz de esta ecuación ubicada 
en el segmento [1, 2]. En este segmento tenemos 

dE 1 1 

ln ma 


Por lo tanto, podemos emplear el método de aproximaciones sucesivas, 
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Al hacer x, = l, tenemos x4 =~ p(x«) = 1,225. Por consiguiente, con un 
error inferior a 0,001 tenemos x = 1,225, 

7. La ecuación dada tiene la raíz x=0. De la figura 26 se ve que 
la segunda raiz de la ecuación es positiva. Por esto, ella satisface la 


ecuación x= |enx. Aqui, p(x) = 


Ya que 
1 cF 1 
o(3)- senz = (04794 > 5 
4 011) = Sen 1 = VÕIS < 1, 


á 1 
entonces, la ecuación tiene una raiz en el segmento [5 5] En este 


segmento tenemos 


t 
cos5 


lg'ixi < 
2 jsn 


y por razón el método de aproximaciones sucesivas converge, Al hacer 
xy = 1, tenemos x7 = Q(x,) > 0,8768, Quiere decir que con un error inferior 
a 0,0001 la segunda raíz de la ecuación es igual a 0,8768. 

8. Se resuelve igual que el anterior. Escribamos la ecuación en forma 


0,8703 
1,3846 


EL 


3 
x= Ysenx, 


ul 


hacemos x, = 1 y obtenemos xe = p (xe) = 0,9286. Por esto, con un error 
inferior a 0,001 una de las raíces de la ecuación es igual a 0,9286. 
Como ambos miembros de la ecuación son las funciones impares, tenemos 
una raíz más, a saber, 09286, La tercera raiz es 0, 


A a x+t 
9. La ecuación podemos escribirla en la forma q =senx 


Fig. 27 


—1/2 < x < 1/2, De la fig. 27 se deduce que esta ecuación tiene la única 
raíz que se encuentra entre O y 2/2. En este caso 


tx) = arcsen S s 


gix s 
yiz +1 


En el segmento [e 3] tenemos |9'()1 < ae 1. Al hacer xı = 0, halla- 
Y 

mos xa = (xa) = 0,3422, Por lo tanto con un error inferior a 0,0001 

tenemos x = 0,3422. 

10, Puesto que cos0 = 1, cosi >0, entonces la ecuación x = cosx 
tiene la raiz en el segmento [0, 1]. Ya que |ọ'(x)| < sen 1 < 1. se puede 
emplear el método de aproximaciones sucesivas. Al hacer x, = 1, obtenemos 
que con un error inferior a 0,0001 x = 0,7391. 

11. La fig. 28 muestra que las raíces positivas de la ecuación son 
próximas a los puntos de intersección del gráfico de la función y= cosx 
con el eje de abscisas y se ubican a la derecha de Jos puntos de intersección 


del tipo 3+ (2k + 1)x, y a la izquierda de los puntos de intersección del 


tipo 3 + 2km. Para hallar la solución ubicada cerca del punto x ajta, 
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n i P 
hagamos x— nī =3=y La ecuación tomará la forma 


n 1 PE 
A r a 
n seny 
cos(y+nn +3 


Fig. 28 


kd r " 
En este caso, 3 sys 3: Y. Por lo tanto, la ecuación puede escribirse 


en la forma 
y =(=1)"*' arcsen ——— . 
y+nn+ pa 
i 2 
Aqui, 
i 
PL) =(— 1)" arcsen x 
v+nr+5 
y y 


a E, e 
7 7 

o rnrts (2900 +5) 

Está claro que cerca del punto y =0 tenemos |y'(y)| <q < l, y por esta 


razón podemos emplear el método de aproximaciones sucesivas. Hallemos 
la solución para n = 1 con un error inferior a 0.001. Sea yy =0 Entonces, 


y2 = píy2) = 0,204. Por consiguiente, y = 0,204; x =} +y 4917 
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Para hallar la primera raiz negativa, hagamos n= — 1 Obtendremos 
la ecuación 


y = arcsen - 


3 
Hagamos yo = 0. Entonces 
Yo = P (y6) = —0,503. 
De este modo, y = —0,503, y, por lo tanto, x = —2,074. 


Con mayores valores de |n| el método de aproximaciones sucesivas 
da la fórmula aproximada para y: 


y = =p ak ala 
y = Qlyo) = 1—17 *' arcsen z a+ i 
nt+5 


Por esta razón 


E Al 
xazin + 1) + TEO 


12. Al hacer x, =0, encontramos que x, = p(x3) = 1,088. Por consi- 
guiente, con un error inferior a 0,00) x = 1,088 
13. Primero resolvamos la ecuación x = |, Ig(x + 2). Tenemos p (x) = 
= gix +2), y por esto 
s lge 
gids —— 
2(x +2) plglx +2) 


Como (0) = } Íg2 > 0, p(1)= 183 < 1, entonces la ecuacion tiene 
una raíz en el segmento [O, 1]. En este segmento se cumple la desigual- 
dad |¢'(x)] <q <1. Por lo tanto, puede emplearse el metodo de apro- 
ximaciones sucesivas. Al hacer x, =l, tenemos xs = (xs) = 0,6507, 


Quiere decir que la raíz de la ecuación x = | Íglx + 2) con un error infe- 
rior a 0,0001 es igual a x = 0,6507 
Examinemos la ecuación 


x=- jigi +2). 


Aquí, 
p= — jigix +2) 


Puesto que p(0)= — | Íg2= —0.55. s(- Y) =- 1Íg13 =—042. dicha 


ecuación tiene la raiz en el segmento [— 1,2, 0]. Al hacer x,=0, 
hallamos que xs = pixs} = — 0,4397. Asi pues, con un error inferior a 
0,0001 tenemos x = — 0.4397. 


14. Una de las raices de la ecuación es v =0. Para hallar la segunda 
raiz, escribamos la ecuación en la forma x= + | intx + 1). Para la ecua- 
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ción x = Yln(x + 1) tenemos 


o6) = [m3 sE, 
e) =/m2<! 


Esto significa que la ecuación tiene una raiz en el segmento (1/2, 1]. 


1 
Puesto que g'(x) = ————=—. entonces en el segmento [1/2, 1 
2ix + 1) n + D) 


tenemos |ọ'(x)] <4 <1. Hagamos x, =1, entonces xg = P(x) = 0,7469. 
Asi pues, con un error inferior a 0,0001 tenemos x = 0,7469. La ecuación 
x= — nix +1) no tiene raices, a excepción de x=0. Asi, x=0 ó 
x = 0,7469. 


15, Escribamos la ecuación en la forma 


x= y4=Inx 
Aquí, 1 
pl) =2, 9) = /4-1n2, g'(x) = — 
y 2xJ/4=Inx 
4- 


Puesto que p(1)=2, p(2) = 4 — In2 <2, entonces la ecuación tiene una 
raíz en el segmento Ñi, 2). La fig. 29 muestra que la ecuación no 


Fig. 29 


tiene otras raíces. Al hacer x, =2, obtenemos 
xa = 0 (xa) = 1,841. 
Asi pues, con un error inferior a 0,001 tenemos x = 1,841. 
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16. Escribamos la ecuación en la forma 


x=2-inx. 
Aquí. etx)=2=Inx. 9'00==2, De la fg 30 se desprende que la raíz 


de la ecuación se ubica en el segmento [1, 2]. En este segmento | y'(x)| < 
<i. Al hacer x, =1.5. encontramos Xi; =Q(x,3) = 1,557. Asi pues, 
x = 1,557 con un error inferior a 0,00t. 

17, De la fig. 31 se desprende que la ecuación sólo tiene una raiz 
negativa. Escribamos la ecuación en la forma 


xa- pet. 


Fig. 30 Fig. 31 


Entonces, 


p)=- PAZ, — 


: 244+2 


pl==-=|2+e7!= 1,54; p(-2)=-]24+0e?x 1,46, 


Por lo tanto, la raiz se ubica en el segmento [-1, —2]. Al hacer 
xı = — 1, encontramos que x, = p(x¿) = — 1.492. Por consiguiente, con un 
error inferior a 0.001 tenemos x= = 1,492. 


18. Está claro que una de las raices de la ecuación es x, = 10. 
Para hallar la segunda raíz. escribamos la ecuación en la forma x = 10%! 
Aquí. 91) = 10, g'(x)=0.1-10%*In10. Además, (1) =10% >1, 
(2) = 10%? < 2. Por esto. la ecuación tiene la raíz en el segmento [i, 2]. 
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¡q_ > A E 02 e —KÉ— 


En este segmento 
L9'(x)| < 0,1- 10%? -In 10 = 0,37 < 1, 


Ahora podemos aplicar el método de aproximaciones sucesivas. Al hacer 
xy =2, encontramos que xy = p(x)) = 1,372. Asi pues, con un error inferior 
a 0,001 tenemos x = 1,372. 

19. De la fig. 32 se deduce que la ecuación tiene raíces en cada uno 
de los segmentos 


[i+ Hn+ ox), md... 


2 


con la particularidad de que estas raices se ubican en las mitades derechas 
de los segmentos. Para hallar la primera raíz positiva, hagamos la susti- 


tución x= m y. La ecuación tomará la forma 
3r 
aail) 


de donde, ya que 0 < y < x, tenemos 


y= arae[ie(32-s)}. 
o0) = wcefis(3x -»)] 


Aqua, 
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En el segmento [0, n} hay raiz de nuestra ecuación, y en este segmento 
|[@'0)] < 1. Hagamos uso del método de aproximaciones sucesivas. Al hacer 
Yı =0, tenemos y4 = plya) = 1,059. Esto significa que con un error infe- 
rior a 0,001 tenemos y = 1,059, y, por lo tanto, x = 3,654, 


Para hallar la segunda raíz positiva, hacemos x = A n — y. La ecuación 


toma la forma 6 
- y= aca 15(5% - »)] 


Al hacer y, =0, obtenemos y4 = ¢ (y4) = 0,876. Por lo tanto, con un error 
inferior a 0,001 tenemos y = 0,870 y x = 6,984. 
20. De la fig. 33 se deduce que la ecuación tiene una raiz ubicada 
L 


entre O y 1. Tenemos g(x) == ye. 010) =—<ye"* En el segmento 


Pig 33 


[0. 1] se cumple la desigualdad | 91 < > que hace posible el empleo 


del método de aproximaciones sucesivas. Al hacer y, =0, hallamos 


Xa S Q (xa) 70,095 Por lo tanto, con un error inferior a 0,0001 tenemos 
x = 0091. 
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21. Sea 
NM0=-5+1 
Entonces 
$) = 3x = 5, f(x) = 6x. 


Según la fórmula de Newton, tenemos 


2—5B,+1 
bmp PEL 


Componemos la tabla de los valores de la función 


x -2 -1 g 3 
fix) =I 3 1 13 


Esta tabla muestra que la ecuación x? — Sx + 1 = 0 tiene las raices en los 
segmentos [-3, -2], [0, 1), [2, 3]. 

Hallemos primero la raíz ubicada en el segmento [ 3, —2]. Puesto 
que en este segmento f”(x) <0, entonces, como primer valor tomamos 
Po = —3 (ya que f(Bo) = —11 es un número negativo). Tenemos 
(EN 

3-3-5 > 


B =-3- -=25. 


Continuando los cálculos, encontramos que B, = Ba = —2,331. Resulta que 
con un error inferior a 0,001 la raíz de la ecuación en el segmento [—3, —2] 
es igual a — 2,331. 

Ahora hallemos la raiz que se encuentra en el segmento [0, 1]. 
Aquí tenemos f”(x)>0. Por ello, hacemos fo =0. De aquí obtenemos 


0-5:0+1 


02, Bs = Ba = 0202. 


B=0- 


Con un error interior a 0,001 tenemos x = 0,202. 
Por fin, para hallar la raiz en el segmento [2,3] hacemos Bo =3 
y tenemos å 
P-SI41 
=3= 200: 
Es 3-5 


Continuando los cálculos, encontramos que Ba = Ps = 2.128. Por esto, con un 
error inferior a 0,001 la raíz es igual a 2,128. Hemos, pues, encontrado 
tres raíces: x; = —2,331, xz = 0,202; x, = 2,128. 
Resolvamos esta misma ecuación por el método de cuerdas perfec- 
cionado. En el segmento [-3, —2] encontramos 
SAN OON E 
A E 
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Como en este segmento f”(x) <0, la curva gira su convexidad hacia la 
parte positiva del eje de las y, por lo cual a, se determina según la 
fórmula 


a= If 


A continuación hallamos 
2 a E 23. 


Esto coincide, con un error inferior a 0,001, con el valor de x hallado 
antes. 
Del modo igual resolvemos la ecuación en los segmentos [0,1] y 
fa, 31 
22. Aquí tenemos 
fix) =x -9x? + 20x— 11, 


$ (x)=3x? — 18x +20, 
f(x) =6x ~ 18 =6(x —3). 


Componemos la tabla de la función f(x); 


y Ó [2 3 ECTS TS 
fo | -n A 3 -5 -u| -u| 1 


Las raices de la ecuación están en los segmentos [0, 1], [2, 3], [5. 6] 

En el segmento [0, 1] hacemos Po = 0 y hallamos que B, = fs = 0,834. 
En el segmento (2, 3] hacemos Bo = 3 y encontramos que Bz = B, = 2,216. 
En el segmento (5, 6] hacemos B, = 0 y encontramos que Ba = [Ps = 5,249 
Encontramos (con un error inferior a 0,001) tres raices de la ecuación 


xı =0,834; xz = 2,216; x, = 5,249 


23. Aquí, f(x) = x? — 3x? = 3x + 11, f (x)= 3x7 — 6x -= 3, J") = 
= 6x — 6 = (x — 1). Componemos la tabla de los valores para f(x): 


y 2 
Its) [to 1 


La ecuación tiene una raíz real ubicada en el segmento [-2, —1]. Para 
hallar esta raiz, hacemos fo = —2. Obtenemos que Bz =f, = — 1,847 
Por esto, con un error inferior a 0,001 tenemos x = — 1,847. 
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24. Aquí f (x}= x$ + 5x + 1, f'(x)= 5x* + 5, f"(x)=20x7. La tabla 
de los valores de la función f(x): 


Por esto, la ecuación tiene una raíz en el segmento [—1, 0].. Hacemos 
Bo = — 1. Tenemos, con un error inferior a 0,0001, fi, =P, = —0,1999, 
Por eso, con la precisión indicada tenemos x = —0,1999. 


25, Aquí, f(x) = sen x + x—1,f'(x) = cosx + 1, f” = — sen x. La tabla 
de los valores de la función: 


Fotr anh] 


La raíz se encuentra en el segmento [0, 1}. Hacemos Bọ =0 y hallamos 
que Pz = B, = 0,5110. Por esto, con un error inferior a 0,0001, tenemos 
x = 0,5110. 


10 
í, =x? - = - —— "e (x) 
26. Aqui, f(x) = x? — 10lg(x) — 3, o) = 2x anio” f'a + 


-r2 La tabla de los valores de la función f(x): 
x*In10 


AAA AAA 
Ga a A) 


Las raices de la ecuación se encuentran en los segmentos [0,5, 1] y [2, 3]. 
En el segmento (0,5, 1] hacemos Bo = 0,5 y obtenemos que Bı = Py = 0,535. 
Por esto, la raiz correspondiente, con un error inferior a 0,001, es igual a 
0,535. En el segmento (2, 3] hacemos Pp = 3 y hallamos que Bz = Py = 2,705. 

La ecuación tiene dos raíces: x, = 0,535; x, = 2,705. 

27. En el sistema a) hacemos xo = 0, yo = Ô, zo = 0 y después de unas 
cuantas aproximaciones encontramos que con un error inferior a 0,001 
tenemos x = 1,021, y = 2,150, z = 3,072 

En el sistema b) hacemos xo =0, yo = 0, y después de unas cuantas 
aproximaciones obtenemos (con un error inferior a 0,001) que x = 0,520, 
y = 0,310. 

En el sistema c) hacemos xọ=0, yọ=0 y hallamos x= 1,000, 
y = 2,000. 
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